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Heat  capacity  and  electrical  resistivity  of  graphite  have  been  measured  in  the 
range  2500  to  3600  K.  Thermal  expansion  of  tantalum  has  been  measured  in  the 
range  of  1300  to  3200  K.  Radiance  temperature  (at  653  nm)  of  tungsten  at  its 
melting  point  has  been  determined.  The  enthalpy  of  graphite  has  been  determined 
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Abstract 


Heat  capacity  and  electrical  resistivity  of  graphite  has  been 
measured  in  the  range  2500  to  3600  K.  Thermal  expansion  of  tantalum 
has  been  measured  in  the  range  of  1500  to  3200  K.  Radiance  temperature 
(at  653  run)  of  tungsten  at  its  melting  point  has  been  determined.  The 
enthalpy  of  graphite  has  been  determined  from  273  to  1173  K.  The  thermo¬ 
dynamic  properties  of  diatomic  molecules  at  elevated  temperatures  has  been 
investigated  and  a  formalism  developed  to  include  contributions  from  both 
the  continuum  states  and  the  metastable  states  associated  with  rotational 
barriers.  The  determination  of  barriers  to  internal  rotation  for  inorganic 
species  has  been  reviewed.  Some  studies  of  reactions  of  M(g)  with  S02(g) 
using  infrared  matrix  isolation  techniques  are  described. 
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RESEARCH  ON  THERMOPHYSICAL  PROPERTIES  BY  DYNAMIC  TECHNIQUES 
A.  Cezairllyan,  A.P.  Miiller,  and  M.S.  Morse 

The  progress  in  research  on  thermophysical  properties  by  dynamic 
techniques  has  been  in  the  following  three  areas: 

1.  Measurements  of  the  heat  capacity  and  the  electrical  resistivity 
of  graphite  in  the  temperature  range  2500  to  3600  K. 

2.  Measurement  of  the  thermal  expansion  of  tantalum  in  the  range 
1500  to  3200  K  by  a  transient  (subsecond)  interferometric 
technique. 

3.  Measurement  of  the  radiance  temperature  (at  653  nm)  of  tungsten 
at  its  melting  point. 
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1.  Measurements  of  the  Heat  Capacity  and  the  Electrical  Resistivity 

of  Graphite  in  the  Temperature  Range  2500  to  3600  K 

Disagreements  exist  between  the  limited  number  of  experimental 
results  for  the  thermal  properties  of  graphite  above  about  1500  K  reported 
in  the  literature,  ^tost  of  the  reported  measurements  were  performed  on 
different  grades  of  graphite,  which  complicates  their  evaluation  and 
poses  the  question  whether  the  differences  were  due  to  measurement 
errors  or  were  really  Indicative  of  differences  in  the  graphite  grades. 

As  an  attempt  to  elucidate  this  problem,  a  program  was  initiated 
for  the  systematic  and  accurate  measurement  of  selected  properties  of 
various  grades  of  graphite  at  high  temperatures.  The  results  of  the 
first  phase,  namely  the  measurements  of  heat  capacity,  electrical 
resistivity  and  hemispherical  total  emittance  of  three  grades  of  graphite 
in  the  temperature  range  1500  to  3000  K  and  in  vacuum  were  presented  in 
an  earlier  publication  [1].  The  results  on  heat  capacity  were  within 
2  percent  for  all  three  grades  of  graphite  (Poco  AXM-5Q,  Poco  DFP-2,  and 
pyrolytic)  in  the  above  temperature  range.  The  high  vapor  pressure  of 
graphite  had  limited  the  performance  of  measurements  to  temperatures 
below  3000  K. 

The  objective  of  the  present  investigation  is  to  extend  the  measure¬ 
ments  to  temperatures  above  3000  K.  In  order  to  achieve  this,  experiments 
were  conducted  with  the  specimen  in  a  pressurized  gas  (argon)  environment 
at  4  MPa  (about  40  atm.).  The  details  of  the  dynamic  technique  used  for 
the  measurements  are  given  in  the  literature  [2,3].  The  major  addition 
in  the  present  work  is  the  high  pressure  cell  which  has  feed-throughs 
for  proper  electrical  connections  to  the  specimen  and  two  optical  windows 
for  pyrometer  sighting  and  for  optical  alignment.  A  schematic  diagram 
of  the  pressure  cell  and  the  specinren  support  stand  is  shown  in  Fig. 

1.1. 

The  measurements  were  performed  on  5  tubular  specimens  (3  of 
Poco  DFP-2  and  2  of  Poco  AXF-9Q1)  fabricated  from  graphite  rods  which, 
according  to  the  manufacturer,  were  99.99+  pure.  Three  experiments  were 
conducted  on  each  specimen.  The  first  experiment  was  performed  with  the 


specimen  in  an  argon  atmosphere  at  atmospheric  pressure  in  the  temperature 
range  2500  to  3000  K.  The  second  experiment,  performed  to  determine  the 
effect  of  pressure  on  the  measurements,  was  conducted  with  the  specimen 
in  an  argon  environment  at  4  MPa  in  the  same  temperature  range.  The 
third  experiment  was  performed  with  the  specimen  in  argon  at  4  MPa  in 
the  temperature  range  3100  to  3600  K.  Typical  heating  rates  for  the 
specimens  were:  1000  K*s~^  at  2500  K,  3000  K*s"^  at  3200  K,  and 
5800  K*s"^  at  3600  K.  Duration  of  the  current  pulses  was  approximately 
270  ms  in  the  lower  temperature  range  and  approximately  300  ms  in  the 
higher  range. 

The  differences  in  the  heat  capacity  results  among  the  two  types  of 
graphite  were  small  and  the  data  from  all  15  experiments  (on  five  specimens) 
were  combined  and  fit  to  a  polynomial  function  in  temperature.  The  best 
fit  was  determined  by  the  standard  deviation  was  obtained  with  a  linear 
function.  The  function  for  heat  capacity  (standard  deviation:  0.7%) 
that  represents  the  results  in  the  temperature  range  2500  K  to  3600  K 
is: 

Cp  =  22.46  +  1.371  X  10"^  T  (1) 

where  T  is  in  K  and  Cp  in  J*mol"^*K”^.  The  heat  capacity  of  graphite 
computed  using  Eq.  (1)  is  given  in  Table  1.  In  the  computations  of  heat 
capacity  the  atomic  weight  of  graphite  was  taken  as  12.011. 

Electrical  resistivity  of  the  graphite  specimens  was  determined 
from  the  same  experiments  that  were  used  to  calculate  heat  capacity. 

However,  unlike  heat  capacity,  there  were  differences  in  the  electrical 
resistivity  of  the  two  grades  of  graphite.  The  results  for  each  grade 
were  separately  fit  to  a  linear  function  in  temperature.  The  functions 
for  electrical  resistivity  that  represents  the  results  for  each  grade  of 
Poco  graphite  in  the  temperature  range  250G  to  3600  K  are: 

Poco  (DFP-2)  graphite  (standard  deviation:  0.3%) 

p  =  618.10  +  0.17349  T  (2) 

Poco  (AXF-9Q1)  graphite  (standard  deviation:  0.4%) 

p  =  618.50  +  0.15946  T  (3) 

where  T  is  in  K  and  p  is  in  pn*cm.  The  electrical  resistivity  of  each 
grade  of  graphite  computed  using  Eqs.  (2)  and  (3)  are  given  in  Table  1. 
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i  Table  1 

! 

i  Heat  capacity  and  electrical  resistivity  of  graphite 


Temperature 

(K) 

Heat  Capacity 
(J*moT^-K‘\) 

Electrical  Resistivity 
(yJi'Cm) 

AXF-9Q1 

2500 

25.89 

1051.9 

1017.2 

2600 

26.02 

1069.2 

1033.1 

2700 

26.16 

1086.6 

1049.0 

2800 

26.30 

1103.9 

1065.0 

2900 

26.44 

1121.3 

1080.9 

3000 

26.57 

1138.6 

1096.9 

3100 

26.71 

1156.0 

1112.8 

3200 

26.85 

1173.3 

1128.8 

3300 

26.98 

1190.7 

1144.7 

3400 

27.12 

1208.0 

1160.7 

3500 

27.26 

1225.4 

1176.6 

3600 

27.40 

1242.7 

1192.6 
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Estimates  of  errors  In  measured  and  computed  quantities  lead  to  the 
following  estimates  of  possible  maximum  errors  in  the  properties:  heat 
capacity,  1%  at  2500  K  and  increasing  to  Z%  at  3600  K;  electrical 
resistivity,  1%  at  2500  K  and  increasing  to  2%  at  3600  K. 

Only  a  limited  number  of  experimental  heat  capacity  results  for 
graphite  at  temperatures  above  2000  K  were  found  in  the  literature.  In 
addition  to  this,  most  of  them  were  for  different  grades  of  graphite. 

Thus,  it  is  questionable  whether  there  is  a  common  base  for  comparison. 

In  order  to  provide  continuity,  results  reported  in  the  literature  at 
temperatures  above  2000  K  are  presented  in  Fig.  1.2. 

The  heat  capacity  results  of  the  present  work  are  in  a  reasonably 
good  agreement  (about  1%  on  the  average)  with  those  reported  earlier  [1] 
for  the  same  grade  of  graphite  in  the  overlapping  temperature  range  2500 
to  3000  K).  The  smoothed  results  of  Rasor  and  McClelland  [4]  are 
average  values,  as  given  by  the  authors,  representing  data  on  four 
grades  of  graphite  (3474  D,  7087,  6BH,  and  GBE).  Their  actual  data  have 
a  scatter  of  approximately  5  percent  about  their  smooth  curve.  In  the 
temperature  range  2500  to  3500  K  the  results  of  Rasor  and  McClelland  are 
on  the  average  about  3  percent  higher  than  those  of  the  present  work. 
However,  at  temperatures  above  3500  K  the  values  reported  by  Rasor  and 
McClelland  increase  rapidly.  The  results  reported  by  Sheindlin  et  al. 

[5]  are  in  agreement,  within  1%  on  the  average,  with  those  reported 
earlier  [1]  in  the  temperature  range  2000  to  2500  K.  However,  above 
2500  K,  their  results  indicate  a  sharp  increase  in  heat  capacity  reaching 
a  value  of  29.6  J*mol"^*K”^  at  3000  K  which  is  about  11  percent  higher 
that  the  present  work  result  and  a  value  of  42.7  J*mol'^*K"^  at  3600  K 
which  is  about  56  percent  higher  than  the  present  work  result.  The 
sharp  increase  in  the  heat  capacity  reported  in  the  above  two  investigations 
is  in  contrast  to  the  results  of  the  present  work.  There  is  no  obvious 
explanation  for  this  difference,  except  that  graphite  evaporation  might 
have  affected  the  results  reported  in  the  literature. 

Electrical  resistivity  of  graphite  is  sensitive  to  the  graphite 
grade  and  also  shows  considerable  variations  for  specimens  from  different 
lots  of  the  same  grade.  Thus,  the  only  meaningful  comparison  of  the 
present  work  results  is  with  those  measured  earlier  on  the  same  grade  of 
graphite  [1].  It  may  be  seen  from  Fig.  1.3  that  the  agreement  is  better 
than  2%  in  the  overlapping  temperature  ran  (2500  3000  K). 


Figure  1.3.  Electrical  Resistivity  of  Graphite  at  High  Temperatures 
Reported  in  the  Literature. 
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With  the  present  system  and  operating  conditions  It  was  not  possible 
to  extend  the  measurements  to  temperatures  significantly  above  3600  K. 
Above  3600  K,  the  smooth  heating  curve  (temperature  versus  time),  as 
obtained  from  the  pyrometer  output,  showed  discontinuities  typical  of 
Interference  of  graphite  vapor  with  radiation  from  the  specimen.  In 
order  to  extend  the  measurements  to  temperatures  above  3600  K,  work  Is 
planned  to  perform  the  measurements  with  the  specimen  under  an  inert  gas 
environment  at  pressures  higher  than  that  used  so  far,  at  least  up  to 
20  MPa. 


References 

1.  Cezalrllyan,  A.,  and  Righini,  F. ,  Rev.  Int.  Htes.  Temp,  et  Refract., 
J2,  124  (1975). 

2.  Cezairliyan,  A.,  J.  Res.  Nat.  Bur.  Stand.  (U.S.)  75C,  7  (1971). 

3.  Cezairliyan,  A.,  Morse,  M.S.,  Berman,  H.A.,  and  Beckett,  C.W., 

J.  Res.  Nat.  Bur.  Stand.  (U.S.)  7^,  65  (1970). 

4.  Rasor,  N.S. ,  and  McClelland,  J.D. ,  J.  Phys.  Chem.  Solids  1^,  17 
(1960). 

5.  Sheindlin,  A.E.,  Belevich,  I.S.,  and  Kozhevnikov,  I.G.,  High 
Temperature  1£,  897  (1972). 


There  has  been  an  increasing  need  in  recent  years  for  thermal 
expansion  data  on  refractory  materials  at  temperatures  above  the  limit 
of  accurate  steady-state  experiments.  Steady-state  techniques  such  as 
push-rod  dilatometry,  x-ray  diffractometry  or  twin- telemicroscope 
methods  are  generally  limited  due  to  problems  associated  with  heat 
losses,  chemical  reactions,  evaporation,  etc.,  which  become  especially 
severe  at  temperatures  above  2000  K.  To  overcome  these  difficulties,  a 
unique  transient  interferometric  technique  has  been  developed  recently 
[1,2]  at  the  National  Bureau  of  Standards  for  the  purpose  of  measuring 
thermal  expansion  of  solids  at  temperatures  primarily  in  the  range 
1500  K  and  the  melting  point  of  the  specimen. 

The  basic  method  involves  rapidly  heating  the  specimen  from  room 
temperature  to  the  maximum  temperature  of  interest  in  less  than  one 
second  by  the  passage  of  an  electrical  current  pulse  through  it,  and 
simultaneously  measuring  the  specimen  temperature  by  means  of  a  high¬ 
speed  photoelectric  pyrometer  and  the  shift  in  the  fringe  pattern 
produced  by  a  Michel  son-type  interferometer.  The  polarized  beam  from  a 
He-Ne  laser  in  the  interferometer  is  split  into  two  component  beams,  one 
which  undergoes  successive  reflections  from  optical  flats  on  opposite 
sides  of  the  specimen,  and  one  which  serves  as  the  reference  beam. 

In  this  report,  we  briefly  summarize  the  results  of  measurements 
performed  during  FY81  on  the  thermal  expansion  of  tantalum  at  temperatures 
between  1500  and  3200  K.  The  purpose  of  these  experiments  was  twofold: 

(1)  to  assess  the  operational  characteristics  of  the  interferometric 
system  over  a  large  temperature  range,  particularly  at  temperatures 
approaching  the  melting  point,  and  (2)  to  obtain  accurate  expansion 
values  for  tantalum  at  temperatures  where  considerable  disagreement 
exists  among  data  reported  in  the  literature. 

The  experiments  were  performed  on  four  specimens,  each  fabricated 
in  the  form  of  a  precision-machined  tube  (with  optical  flats)  containing 
a  small  sighting  hole  for  direct  pyrometric  measurement  of  blackbody 
temperature.  In  a  typical  experiment,  the  specimen  was  resistively 
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pulse-heated  in  a  vacuum  environment  (^-l  mPa)  from  room  temperature  to 
the  desired  temperature  in  about  800  ms.  The  analog  signals  from  the 
pyrometer  and  the  interferometer  were  recorded  during  the  experiment  by 
means  of  digital  storage  oscilloscopes.  The  recorded  data  was  then 
transferred  to  a  minicomputer  for  subsequent  computation  of  specimen 
temperatures  and  fringe  shifts. 

The  linear  expansion  of  a  specimen  was  determined  at  each  recorded 
temperature  from  the  cumulative  fringe  shift  which,  for  the  range  1500 
to  3200  K,  varied  from  about  165  to  480  fringes.  The  expans ion/temperature 
data  pairs  for  the  four  specimens  were  combined  and  then  fitted  by  a 
polynomial  function  of  temperature  by  means  of  the  least-squares  method. 

The  function  that  represents  the  results  for  linear  thermal  expansion  of 
tjit^talum  (standard  deviation  0.2%)  in  the  temperature  range  1500  to 
3200  K  is 

=  5.141  X  lO-**  +  1.445  X  lO-^  T  +  4.160  x  lO'^  T^ 

-  1.309  X  10"^^  T^  +  1.901  X  10'^®  T**  (1) 

where  T  is  in  K  and  is  the  specimen  length  at  20'’C  (nominally  6.1  inn 
for  all  specimens).  The  smoothed  results,  as  defined  by  Eq.  (1),  are 
given  (in  percent)  at  intervals  of  100  K  in  Table  1. 

The  results  for  each  specimen  were  also  fitted  by  a  quartic  polynomial 
in  temperature  by  means  of  the  least-squares  method.  As  shown  in  Fig.  2.1, 
the  deviation  of  the  smoothed  results  for  the  individual  specimens  from 
those  represented  by  Eq.  (1)  varies  from  about  0.2%  at  2000  K  to  about 
0.4%  at  3000  K. 

The  major  source  of  error  in  our  results  arises  in  the  determination 
of  specimen  temperature.  It  is  estimated  that  the  error  in  temperature 
is  about  5  K  at  2000  K  and  about  10  K  at  3000  K  which  correspond  to 
uncertainties  in  of  about  0.3%  and  0.5%  at  the  respective  temperatures. 

The  maximum  uncertainty  in  the  fringe  count  is  believed  to  vary  from 
about  0.5  fringe  at  2000  K  to  about  1  fringe  at  3000  K  which  correspond 
to  an  uncertainty  in  of  approximately  0.2%  over  the  entire  temperature 

range.  A  further  uncertainty  of  about  0.2%  arises  in  determining  the 
specimen  length  at  "room"  temperature  (20°C).  Therefore,  the  maximum 
error  in  our  reported  values  of  linear  expansion  is  estimated  to  be 
about  1%  at  2000  K  and  not  more  than  2%  at  3000  K. 
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Table  1.  Smoothed 

results  for  linear  thermal  expansion 

of  tantalum. 

Temperature, 

At/tg. 

Temperature, 

Al/tp, 

K 

% 

K 

% 

0.859 

2400 

1.616 

0.936 

2500 

1.710 

1700 

1.015 

2600 

1.807 

1800 

1.096 

2700 

1.908 

1900 

1.178 

2800 

2.012 

1.261 

2900 

2.121 

1.347 

3000 

2.234 

1.434 

3100 

2.353 

1.524 

3200 

2.478 

Figure  2.1.  Deviation  of  the  smooth  results  for  the  individual  specimens 
from  the  least-squares  fit  [represented  by  Eq.  (1)]  to  the 
combined  data  for  the  four  specimens. 


reported  in  the  literature. 
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The  results  of  the  present  work  are  compared  in  Fig.  2.2  with 
expansion  values  reported  in  the  literature.  Worthing  [3],  Rasor  and 
McClelland  [4]  and  Conway  et  a1.  [5]  obtained  their  expansion  data  using 
the  twin-telemicroscope  method.  The  values  of  expansion  reported  by 
Edwards  et  al.  [6]  and  Waseda  et  al.  [7]  were  measured  by  x-ray  diffraction 
techniques.  Amonenko  et  al.  [8]  obtained  their  data  by  conventional 
dilatometry.  The  data  disagree  by  as  much  as  10%  throughout  the  entire 
temperature  range  but  do  not  show  any  significant  bias  towards  a  given 
steady-state  measurement  technique.  The  data  reported  by  Worthing  and 
by  Edwards  et  al.  are  in  rather  good  agreement  with  the  results  of  the 
present  work.  The  expansion  values  determined  by  Conway  and  Losekamp 
agree  with  those  of  the  present  work  within  the  combined  experimental 
errors  of  the  two  investigations.  However  the  expansion  results  reported 
by  the  remaining  investigators  are  considerably  higher  than  our  results. 

In  summary  the  results  of  the  present  study  show  that  our  transient 
interferometric  technique  is  capable  of  measuring  thermal  expansion  of 
metals  at  temperatures  between  1500  K  and  their  melting  points  with  an 
estimated  error  of  less  than  2%. 
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3.  Measurement  of  the  Radiance  Temperature  (at  653  nm)  of  Tungsten 

at  its  Melting  Point 

With  the  increased  interest  in  high  temperature  applications, 
accurate  temperature  measurements  in  high  temperature  systems  has  become 
important.  A  major  problem  is  the  difficulty  associated  with  reference 
points  in  temperature  calibrations.  The  best  conventional  reference 
points  are  blackbody  furnaces  operating  at  the  melting  point  of  selected 
metals.  Because  of  the  serious  problems  associated  with  systems  exposed 
to  high  temperatures  for  prolonged  periods  of  time,  the  practical  limit 
for  blackbody  furnaces  is  2000  K.  During  recent  years,  we  have  been 
investigating  the  possibility  of  establishing  secondary  reference  points 
at  high  temperatures  based  on  the  surface  radiance  temperature  of  metals 
during  melting.  Our  earlier  work  on  radiance  temperature  measurements 
at  the  melting  point  of  selected  refractory  metals,  such  as  niobium  [1], 
molybdenum  [2],  tantalum  [3],  indicates  reliability  and  reproducibility 
of  such  measurements. 

The  present  work  Involved  the  measurement  of  radiance  temperature 
(at  653  nm)  of  tungsten  at  its  melting  point.  Since  tungsten  has  the 
highest  melting  point  among  the  elements,  the  results  are  of  particular 
interest  to  the  high  temperature  field. 

The  measurements  were  performed  on  23  tungsten  specimens  of  99.9'’' 
purity.  Before  the  experiments,  the  surfaces  of  a  number  of  specimens 
were  treated  with  an  abrasive.  Two  different  grades  of  abrasive  were 
used  yielding  two  surface  roughnesses  (approximately  0.2  and  0.4  pm  in 
roughness). 

All  the  experiments  were  performed  with  the  specimen  in  an  argon 
environment  at  atmospheric  pressure.  The  heating  rates  for  different 
specimens  were  in  the  range  1100  -  5400  K*s~\  Variation  of  the  radiance 
temperature  as  a  function  of  time  for  two  typical  experiments  is  shown 
in  Fig.  3.1. 
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It  may  be  seen  that,  regardless  of  the  differences  in  the  initial 
surface  conditions,  radiance  temperature  at  the  melting  plateau  is 
approximately  the  same  for  the  two  specimens.  The  same  conclusion  is 
reached  for  all  other  specimens.  For  each  specimen,  the  plateau  temperatures 
were  averaged  to  yield  a  single  temperature  for  that  specimen.  The 
standard  deviation  of  these  averages  ranged  from  0.7  to  0.9  K.  Fig.  3.2 
shows  the  departure  of  the  average  radiance  temperature  for  each  specimen 
from  the  average  for  all  the  specimens.  The  overall  average  for  23 
specimens  yielded  a  value  of  3208  K  for  the  radiance  temperature  (at  653 
nm)  of  tungsten  at  its  melting  point  with  a  standard  deviation  of  0.8  K. 

There  is  no  evidence  of  any  significant  bias  with  respect  to  initial 
surface  roughness  of  the  specimen  or  the  heating  rate  of  the  specimen. 

The  total  error  in  the  measured  radiance  temperature  is  estimated  to  be 
not  more  than  10  K. 

The  results  on  tungsten  have  confirmed  the  earlier  results  on  other 
refractory  metals  that  radiance  temperature  at  the  melting  point  of 
metals  is  reproducible  and  constant  for  a  metal  and  is  virtually  independent 
of  the  initial  surface  conditions  of  the  specimen  and  the  operational 
conditions  (heating  rate,  etc.)  of  the  overall  system.  In  conclusion, 
the  results  suggest  the  possibility  of  using  the  radiance  temperature  at 
the  melting  point  of  selected  metals  for  secondary  calibration  and 
checking  of  optical  temperature  measuring  equipment  in  high  temperature 
systems. 
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Enthalpy  of  POCO  Graphite  AXM-5Q1  From 
273  K  to  1173  K 
by 

David  A.  Ditmars 

1.  Introduction 

Graphite  is  an  important  material  in  high- temperature  technology,  finding 
applications  as  both  a  standard  reference  material  for  high- temperature 
thermophysical  properties  and  as  a  component  in  carbon-carbon  composites 
for  high- temperature  structural  components.  For  this  reason,  accurate 
determination  of  the  thermophysical  properties  of  graphite  has  been  the  subject 
of  several  investigations  in  recent  years  [1-7].  The  most  important  thermophysical 
properties  for  the  purpose  of  heat-transfer  calculations  are  thermal  diffusivity, 
thermal  conductivity,  and  heat  capacity.  These  are  related  to  each  other  in  the 
relation. 


where  k  is  thermal  diffusivity;  k,  thermal  conductivity;  p,  density;  and  Cp, 
specific  heat.  Since  each  of  the  three  thermophysical  properties  contained  in 
this  relation  can  be  measured  independently,  the  possibility  exists  for 
comparing  directly-measured  values  for  any  one  of  these  properties  with 
corresponding  values  calculated,  using  the  above  relation,  from  direct 
measurements  of  the  other  two  properties. 

Graphite  is  obtainable  in  both  natural  and  synthetic  forms.  The  synthetic 
form  is  available  in  many  different  "grades".  These  grade  designations  are 
designations  assigned  by  the  manufacturer  which  specify  in  coded  form  such 
information  on  the  graphite  as  grain  size,  temperature  of  graphitization, 
chemical  purity,  or  other  proprietary  data. 
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The  specific  heat  data  for  all  graphite  grades  has  shown  a  rather 
large  variance,  most  data  lying  in  a  range  +  10%  between  room  temperature 
and  2500  K  [7].  It  has  not  been  made  clear  whether  this  large  variance  is 
directly  related  to  systematic  measurement  error  or  whether  it  indicates  a 
real  dependence  of  thermophysical  properties  on  the  method  of  synthesis. 

For  this  reason,  it  has  seemed  desirable  to  have  a  set  of  data  for  at  least 
one  of  the  thermophysical  properties  cited  above,  extending  over  a  wide 
temperature  range, 

A  very  pure  grade  of  graphite  ("POCO"  graphite)  produced  by  the  POCO 
graphite  company  of  Decatur,  Texas,  and  designated  as  grade  AXM-5Q1  has 
been  chosen  by  NBS  as  a  high-temperature  thermal -conductivity  standard 
reference  material.  Direct  measurement  of  the  heat  capacity  of  this  material 
by  a  pulse-calorimetric  technique  in  the  temperature  range  1500  to  3000  K 
[3,4]  has  been  carried  out  earlier  in  the  NBS  Thermophysics  Division.  Specifi 
heat  data  on  this  same  material  obtained  by  other  investigators  have  also 
been  reported  in  the  recent  literature  [5]  in  the  temperature  range  300  to 
1000  K.  These  data  were  obtained  by  scanning  calorimetry. 

We  have  undertaken  the  measurement  of  the  specific  heat  of  graphite, 
grade  AXM-5Q1  in  the  temperature  range  273  to  1500  K  by  relative-enthalpy 
calorimetry  using  a  precision  Bunsen  ice  calorimeter  and  a  precision 
adiabatic-receiving  calorimeter.  The  aim  is  to  provide  more  accurate 
specific-heat  data  below  1000  K  and  to  bridge  the  present  gap  in  the  specific- 
heat  data  for  this  material  between  1000  and  1500  K.  This  report  gives  the 
preliminary  results  on  enthalpy  data  obtained  in  the  temperature  range  273  to 
1173  K. 
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2.  Samples 

Two  samples  of  graphite  were  chosen  for  enthalpy  measurement.  These 
graphites  are  designated  as  samples  "1"  and  "2".  Sample  1  comprised  two 
specimens  of  AXM-5Q1  graphite,  obtained  from  the  lot  of  material  reserved 
for  the  NBS  thermal -conductivity  SRM.  Sample  2  consisted  of  a  single 
specimen  of  AXM-5Q  graphite  (in  principle,  a  grade  of  slightly  lower  purity) 
chosen  from  the  same  lot  which  was  the  source  of  the  graphite  specimens  which 
were  measured  in  the  prior  AGARD  project  of  cooperative  measurements  on  heat- 
transport  phenomena  of  solids  at  high  temperatures  [1].  It  had  been  previously 
noted  [8J  that  there  was  a  small  (<  4  %)  density  variance  throughout  the 
stock  of  NBS  graphite  (sample  1).  Therefore,  the  two  specimens  of  sample  1 
for  enthalpy  measurement  were  chosen  from  rods  "62-end  2"  and  "104-end  2", 
which  showed  a  representative  density  difference.  These  specimens  are 
hereafter  designated  "lH-62"  and  "lH-104".  The  specimen  of  the  sample  2  for 
enthalpy  measurement  is  designated  "2H".  These  three  specimens  were 
machined  as  right  circular  cylinders  for  encapsulation  and  enthalpy  measurement. 
The  measured  specimen  masses  and  calculated  densities  are  given  in  Table  1. 

Table  1 

Enthalpy -measurement  specimens  of  POCO  graphite 


Calculated 

Specimen 

Specimen 

Grade 

Apparent  density 

Mass 

-3 

g  cm 

g 

(lH-62) 

AXM-bQl 

1.738 

6.98047 

(lH-104) 

AXM-5Q1 

1.706 

6.82563 

(2H) 

AXM-5Q 

1.756 

7.05369 
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Chemical  analyses  on  the  present  POCO  graphite  were  performed  by  two 
independent  analytical  laboratories  on  specimens  chosen  from  sample  1;  one 


of  these  laboratories  also  included  a  specimen  of  sample  2  in  their  analyses. 

Specimens  of  both  samples  were  subjected  to  combustion  analyses  by  the 
Galbraith  Laboratories  of  Knoxville,  Tennessee.  These  analysis  specimens 
were  cut  from  portions  of  the  appropriate  sample  materials  which  were  adjacent 
to  the  enthalpy-measuring  specimens.  The  results  of  these  analyses  are  given 
in  Table  2.  Here,  the  specimen  nomenclature  differs  from  that  of  Table  1 
only  in  that  "A"  is  substituted  for  "H". 


Table  2 

Combustion  Analyses  of  POCO  Graphite 
Calculated  _ Composition  -  wt  % 


Specimen 

Apparent  Density 

C 

H 

N 

ash 

g  cm"^ 

(lA-62) 

1.731 

99.66 

<’'0.05 

0.027 

0.29 

(1A-104) 

1.728 

97.77 

<  0.05 

0.017 

1.94 

{2A) 

— 

99.49 

<  0.05 

0.012 

0.48 

Specimens  of  Sample  1  were  examined  in  the  analytical  laboratories  of  the 
Union  Oil  Company  of  California®.  These  specimens  were  chosen  from  different 
rods  of  Sample  1  designated  rods  "13-end  2"  and  "46-end  2".  The  corresponding 
analysis  specimens  are  designated  (lA-13)  and  (lA-46).‘  Table  3  and  Figures  1 
and  2  present  the  results  of  porosimetric,  qualitative  spectrometric  and 
absorbed  gas  analyses  on  these  specimens.  Mercury  porosimetry  analysis 

®  We  are  indebted  to  Dr.  M.J.  Block  of  the  Union  Oil  Company  for  his  generosity 
and  cooperation  in  making  the  facilities  for  these  analyses  available. 
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(Figures  1  and  2)  of  these  specimens  determined  the  total  open  porosity  (OP') 
to  be  0.08114  and  0.08053  cc/g  for  (lA-13)  and  (lA-46),  respectively.  These 
are  converted  to  percent  total  volume  units  by  the  following  equation 
OP  =  (OP*)  (p)  (100), 

where  OP  is  the  percent  open  porosity  and  p  is  the  calculated  density  of  the 
graphite.  The  theoretical  porosity  is  determined  by  the  equation 

and  is  a  measure  of  the  porosity  expected  when  the  theoretical  density  of 
graphite  (2.26  g/cc)  and  the  observed  density  of  the  graphite  artifact  are 
taken  into  account.  The  open  porosity  is  always  less  than  the  theoretical 
porosity  due  to  the  fact  that  some  fraction  of  the  porosity  is  blind  and 
does  not  communicate  with  the  outside  of  the  piece.  The  difference  between 
the  theoretical  porosity  and  open  porosity  is  the  closed  porosity.  Pore 
diameter  and  surface  area  data  are  also  derived  from  porosimetry  testing 
tnd  are  also  listed  in  Table  3. 

Elemental  analyses  performed  on  the  graphite  specimens  were  limited 
to  hydrogen  analysis  and  qualitative  emission  spectroscopy  (ES)  analysis  of 
the  graphite  ash.  These  results  are  also  presented  in  Table  3.  An  attempt 
was  made  to  analyze  the  gas  adsorbed  on  the  graphite.  A  sample  of  graphite 
was  ground  to  a  fine  powder  and  then  sealed  in  an  evacuated  glass  bomb.  The 
sample  was  heated  to  about  200  °C  and  the  off-gas  was  analyzed  via  mass 
spectroscopy.  The  major  components  of  the  off-gas  were  carbon  dioxide,  carbon 
monoxide,  and  water.  Table  3  supplies  further  details  concerning  the  off-gas 
composition. 
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Table  3 

Porosimetric,  Qualitative  Spectrometric  and 
Adsorbed  Gas  Analyses  of  POCO  Graphite 


Specimen  (lA-13) 


Density,  g/cc  (calculated)  1.766 

Porosimetry  Data 

Pore  Volume,  %  volume 


Theoretical 

21.86 

Open 

14.33 

Closed 

7.53 

Average  Pore  diameter,  p 

1.21 

Surface  Area, 

m'^/g 

0.52 

Elemental  Analyses 

Hydrogen,  Wt 

% 

0.005 

Ash,  Qualitative  E.S. 

10% 

Si 

1-10% 

— 

0.1-1% 

Fe,  Ca, 

0.1% 

Mg,  Sr, 

Adsorbed  Gas  Analysis 

Species 

Mole  % 

Species 

CO2 

42.5 

«2 

H2O 

39.6 

^2^2 

CO® 

13.3 

I-C5 

Air 

2.94 

HpS 

CH4 

0.48 

■-C4 

^3^8 

0.43 

':2"6 

SO2 

0.22 

C4H8 

C2H4 

0.13 

COS 

C,H4 

0.11 

a.  A  maximum  of  3%  could  be  N». 

(lA-46) 

1.705 


24.56 

13.73 

10.83 

1.12 

0.64 

0.022 

Si 

V,  Al,  Pb  Fe,  Al,  Cu,  Pb 

Cu,  Cr  Mg,  Cr,  Cu,  Ni,  Mn,  Sr 


Mole  % 
0.10 
0.06 
0.06 
0.03 
0.03 
0.02 
0.01 
0.01 
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3.  Experimental  Procedure 

The  three  graphite  specimens  were  sealed  with  no  prior  treatment  in 
Pt-lORh  capsules  containing  helium  gas  at  a  pressure  of  approximately  0.2  atm. 
Relative  enthalpy  measurements  (Hj-H2y3  ^5)  were  made  on  each  of  the  three 
encapsulated  specimens  in  the  temperature  range  273-1173  K  with  a  Bunsen  ice 
calorimeter.  The  measuring  apparatus  and  techniques  have  been  described 
earlier  [9-11].  Enthalpy  data  were  obtained  at  nine  temperatures  starting 
at  373  K  and  continuing  at  nominal  100  K  intervals.  Specimen  (lH-62)  was 
measured  two  or  three  times  at  each  temperautre  to  enable  estimation  of  the 
precision  of  the  data  series;  the  other  specimens,  not  at  each  temperature 
and  not  as  frequently.  Empty  capsule  data  obtained  previously  [10]  for 
identical  Pt-lORh  capsules  were  used  in  calculating  the  net  contributions 
of  the  graphite  to  the  gross  observed  heats. 

In  addition  to  the  graphite  enthalpy  measurements,  measurements  were 
made  on  a-Al202  (NBS  Standard  Reference  Material  720)  at  selected  temperatures 
to  verify  the  overall  accuracy  of  the  measuring  technique. 

4.  Results 

The  enthalpy  data  on  the  three  graphite  specimens  is  given  in  Table  4  and 
Figure  3.  The  relative  enthalpy  data  of  Table  4  were  calculated  by  subtracting 
from  the  gross  measured  heat  the  capsule  enthalpy,  taken  from  [10],  at  the 
corresponding  temperature  and  dividing  by  the  appropriate  sample  mass  to 
express  the  results  on  a  unit  mass  basis.  For  the  purpose  of  these  calcu¬ 
lations,  the  specimens  were  assumed  to  be  100  percent  graphite.  The  contribution 
of  the  specimen  capsule  to  the  gross  measured  heats  did  not  exceed  25  percent. 
Enthalpy  measurements  on  a-Al202  made  at  373  K,  273  K,  and  1173  K  agreed  with 
earlier  enthalpy  data  on  this  Standard  Reference  Material  to  better  than 
0.15  percent. 


The  smoothed  form  for  these  enthalpy  data  will  be  influenced  also  by 
existing  heat-capacity  data  on  graphite  from  low- temperature  adiabatic  calori¬ 
metry  and  on  future  measurements  in  this  laboratory  on  these  specimens  in 

t 

the  1173  K  to  1500  K  region.  In  Figure  3,  however,  we  have  compared  the 
present  data  with  enthalpies  derived  from  the  smooth  heat-capacity  data 
presented  in  [7]  for  AXM-5Q1  graphite.  Above  373  K,  there  is  no  significant 
difference  in  the  enthalpy  data  for  the  three  specimens.  The  deviations  of 
the  data  for  specimens  (lH-104)  and  (2H)  at  373  K,  while  coming  in  a  temperature 
range  in  which  the  calorimeter  imprecision  has  been  historically  high,  are 
somewhat  larger  than  anticipated.  It  is  not  yet  clear  whether  the  difference 
arises  from  random  error  or  is  due  to  water  or  other  substances  adsorbed  on 
the  specimens. 
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The  Thermodynamic  Properties  of  Diatomic  Molecules  at  Elevated 
Temperatures:  Role  of  Continuum  and  Metastable  States 


Frederick  H.  Mies  and  Paul  S.  Julienne 


ABSTRACT 


A  complete  quantum  formulation  of  diatomic  partition  functions 
uses  the  energy-variation  of  the  elastic  scattering  phase-shift  to 
represent  the  phase  space  associated  v;ith  the  molecular  continuum 
states.  The  resonance  structure  in  the  phase-shift,  due  to  tunnelling 
through  rotational  barriers,  gives  a  rigorous  interpretation  of  the 
metastable  states  which  lie  behind  the  barrier,  and  we  can  justify  the 
need  to  include  such  states  in  the  evaluation  of  thermodynamic 
properties.  However,  we  also  find  that  it  is  inconsistent  to  merely 
include  the  metastable  phase-space  without  considering  the  remaining 
contributions  from  the  continuum.  If  either  component  is  ignored  in 
treating  the  dimer,  then,  of  necessity,  their  presence  will  appear  as 
the  virial  coefficients  which  cause  non-ideal  behavior  for  the  atomic 
fragments  of  the  molecule.  Both  approaches  are  consistent  and  yield  a 
proper  equation  of  state  for  the  vapor.  We  will  show  quantitative 
agreement  between  the  exact  quantal  results  and  the  approximate 
classical  expressions  for  the  partition  function  at  high  temperatures. 
Since  the  classical  theory  automatically  includes  all  the  effects  of  the 
molecular  continuum,  v/e  suggest  that  this  is  both  the  simplest  and  most 
unambiguious  procedure  for  extending  thermodynamic  tables  to  elevated 
temperatures.  Explicit  calculations  are  presented  for  Li2  and  Na2. 


I.  INTRODUCTION 


At  elevated  temperatures  the  equation  of  state  and  other  observable 
equilibrium  properties  associated  with  the  diatomic  molecule  A2  must 
account  for  the  dissociation  processes, 

^  A(oc)  *  A(oc') .  (0 

The  ensemble  is  divided  into  sub-sets  of  molecules  AgCi)  in  specific 
(1)^ -degenerate  electronic  states  |i>.  Each  state  is  characterized  by  an 
interatomic  potential  V^(R)  which  adiabatically  correlates  with  a  unique 
pair  of  atomic  states,^ 

Vj.Ik.'i  -  E^a)*F^.U).  (z) 

■R  o* 

The  internal  partition  function  associated  with  each  attractive  state 
contains  both  bound  Q^(B)  and  continuum  or  dissociative  Q^(C)  state 
contributions.^’^’^’^ 

Q;.  “  ■*"  Qa!  .  (3) 


At  low  temperatures  is  dominated  by  the  finite  set  of  stable,  bound 
rotational  (ji)  -  vibrational  (v)  eigenstates  with  discrete  energies 
0  v/hich  define  Q4(B), 

1  f  ,V  1 


Q^Cb-)  = 


ir 


-  e 


/AT 


It  is  this  contribution  which  normally  defines  the  molar  thermodynamic 
parameters  we  associate  with  the  stable  molecule  A2.  Hov/ever,  at  high 
temperatures,  or  for  predominately  repulsive  potentials,  we  must  expect 


significant,  if  not  dominant,  contributions  from  the  continuum  of  states 
with  e>0  which  lie  above  the  dissociative  limit. 


We  shall  examine  the  exact  quantum  mechanical  expression  for 

0.(C)  =  7  o(6  e  (5) 

*  A  “  -Fae —  • 


Q 

The  density  of  molecular  continuum  states  is  related  to  the  phase-shift 
Us  (,(£)  for  the  scattering  of  atoms  A(a)  +  A(a')  by  the  electronic- 

Q 

rotational  potential. 


+ 


nrn^  ^ 


where  the  relative  kinetic  energy  e  is  measured  with  respect  to  the  R  =  » 
asymptote  for  v^(R), 

^c'  ~  Vji(R)  -  U)  -B^t  U)  ^ 

R->  oo 

In  particular  we  are  concerned  with  the  role  of  metastable 
rotational -vibrational  levels,  e.  >0,  which  are  predicted  to  exist 

1  f  X  ylTi 

3  4  5 

behind  the  rotational  barriers  ’  ’  in  Eq.  (6),  and  we  shall  segregate 
Q^(C)  into  two  components. 


Q/  Cc)  =  Q^CM)  QiCn,  (g) 


where. 


Q/Cm)  = 

JL  'Tri 

This  finite  set  of  predissociating  levels  is  often  combined  with 
Q^(B)  to  define  the  thermodynamic  properties^®’^^  of  the  dimer  A^,  and 
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there  has  been  some  confusion  as  to  whether  this  procedure  is 

legitimate.  We  shall  show  that  this  approach  can  be  justified  on 

rigorous  quantum  mechanical  grounds  in  tenns  of  the  resultant  resonance 
581213 

structure  ’  *  *  contained  in  the  phase-shift  .(e)-  However,  we 

1  9X 

also  contend  that  the  residual  contributions  contained  in  what  has  been 
termed  the  free-continuum  ’  component  Q^(F)  cannot  be  neglected. 

Although  Q^.(B)  and  Q^.{M)  are  always  positive,  Q^(F)  can  be  of 
either  sign  depending  on  v^  and  T,  and  it  is  inconsistent  to  neglect  the 
free  continuum  when  defining  the  properties  of  We  will  propose  that 
the  most  accurate,  and  by  far  the  simplest  course  is  to  use  the  complete 
classical  theory  for  at  high  temperatures.  The  partition  function  is 
directly  related  to  the  second  virial  coefficient  for  potential  v^, 

,  (10) 


where  is  an  abbreviated  notation  for  the  translational  partition 

function  per  unit  volume  for  particles  of  mass  m, 

=-  /  T-TT  m  \ 


=  ( 


)■ 
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The  simple  classical  expression  for  B^,  ’  ’  *  ’ 


By  =  ZTr- 
^  e> 


n 


contains  all  three  contributions  to  Q.,  i.e., 

i-TS^CF). 


c«o 


-/)  (l^) 
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and  therefore  there  is  no  bias  introduced  into  the  evaluation  of  the 

5 

thermodynamic  variables.  We  have  already  shov/n  that  the  classical 
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theory  yields  quantitative  agreement  with  the  quanta!  theory  once  kT 
exceeds  the  vibrational  spacings. 

Equation  (1)  implies  that  a  total  number  of  atoms  N  confined  in 
some  volume  V  at  temperature  T  will  partition  into  two  portions  and 
which  much  satisfy  the  conservation  of  mass 

N  =  +2 


These  entities  which  we  refer  to  as  'monomers'  and  'dimers'  are  defined 
by  the  particular  regions  of  phase  space,  or  partition  functions  (p.f.)^ 
and  (p.f.)y^»  that  we  choose  to  assign  to  the  particles.  Implicit  in 
all  thermodynamic  tabulations^^*^^ is  the  assumption  that  the 
center-of-mass  motion  of  the  atoms  A  and  the  molecules  Ag  behave  as 
ideal  gases,  such  that  the  equation  of  state  is 

r/hr  =  A/^/v 

=  N/v  —  /v . 


US') 


Given  the  equilibrium  constant 

Waa  __ 


N, 


we  can  express  (15)  as  an  infinite  power  series  in  (N/V), 

3 


p  o  c  g  7 

which  is  comparable  to  the  usual  virial  expansion  '  ’  ’  ’ 


AT 


r  ^ 

k  I 


-<C>  (^)  ^ 


(17) 


(If) 
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v/ith  the  result  that 


<B>  ^  c-j) 

Note  however  that  the  molecular  dimer  contributes  to  all  virial 
coefficients  and  <C>  is  negligible  only  if  N/V)«l.  In  fact  Eq. 
(17)  is  not  convergent  unless  (K  N/V)<;l/8,  and  for  strongly  bound 
molecules  it  is  preferable  to  employ  Eq.  (15)  which  describes  both  low 
and  high  temperature  limits  exactly. 

The  decomposition  of  the  internal  partition  function  into  bound 
(B),  metastable  (M),  and  free-continuum  (F)  components  in  Eqs.  (3)  and 
(8)  leads  to  a  similar  separation  of  the  equilibrium  constant 

and  the  dimer  concentration 

CZl) 

in  Eq.  (16).  Thus  we  can  cast  the  equation-of-state  into  several 
equivalent  forms  ^ 


J1  ^ 

CZZou') 

hr- 

V  V 

Czzy) 

V  ^  V  e  '  V 

iZZ<L) 

all  subject  to  the  constraint  imposed  by  Eq.  (14).  Obviously  by 
modifying  the  activity  of  the  moncn'er  concentration,^  from  ideal  in  Eq. 


(22a)  to  an  imperfect  gas  with  a  second  virial  equal  to  in  Eq. 

(22b)  and  ■•^gq(C)  in  Eq.  (22c),  we  can  obtain  different  interpretations 

of  the  residual  molecular  dimer  concentrations.  This  points  to  the 

basic  problem  inherent  in  the  use  of  tabulated  thermodynamic  data  to 

calculate  equilibrium  properties  at  high  temperatures.  Unless  ^^^(F) 

and  Kgq(M)  are  negligible  we  must  be  cognizant  of  the  specific  regions 

of  phase  space  that  have  been  included  in  the  definition  of  the  species 

A2.  This  supports  the  virtue  of  using  the  unambiguous  expression  in 

Eq.  (12)  to  calculate  K  . 

eq 

The  quantum  theory  of  diatomic  thermochemistry  is  discussed  in 
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section  II  with  specific  application  to  the  alkali  dimers  shown  in 
Fig.  1.  Scattering  theory  is  used  to  analyze  the  metastable  and 
continuum  contributions  to  the  partition  functions  in  section  III, 
Explicit  calculations  of  equilibrium  constants  for  Li2  and  Na2  are 
presented  and  discussed  in  section  IV.  A  brief  summary  and 
recommendations  for  calculating  high  temperature  thermodynamic  variables 
are  contained  in  Section  V. 


-  R  - 


II.  QUANTUM  THEORY  OF  DIATOMIC  THERMOCHEMISTRY 


The  'free'  atoms  A  and  the  diatomic  molecule  A2  in  Eq.  (1)  are 
uniquely  defined  by  their  partition  functions.  However  we  must  fulfill 
two  constraints  if  we  hope  to  convert  the  resultant  functions  into 
observable  physical  properties  of  the  equilibrated  gas,  such  as  the 
equation  of  state^ 

rr  •'•a  ■  -  .  ▼*  •’Ara  >  -  ■  ■  - 


Hr 
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'is  V 


First  v/e  must  insure  that  the  portions  of  phase  space  assigned  to  N^  and 
in  Eq.  (16)  do  not  overlap.  This  condition  is  automatically 

O 

satisfied  by  the  definition  of  the  phase-shift  we  will  employ  in  Eq. 
(5). 


Second,  we  should  be  convinced  that  together  the  two  partition 

functions  encompass  the  entire  phase  space  accessible  to  the  N  particles 

in  Eq.  (14).  The  interactions  between  A+A2  and  A2+A2  must  ultimately 

modify  the  activities  of  N^  and  at  high  pressures.  However  our 

concern  is  with  strongly  bonded  dimers  with  deeply  attractive  ground 

states  as  represented  by  the  alkali  potentials  in  Fig.  1,  and  it 

is  reasonable  to  assume  there  is  a  wide  regime  where  Eq.  (23)  is  an 

accurate  representation  of  the  system.  In  any  case,  our  goal  in  this 

paper  is  not  to  perform  a  complete  analysis  of  the  alkali  vapors,  but 

rather  to  obtain  a  complete  understanding  of  the  entity  we  commonly 

refer  to  as  the  diatomic  molecule  intent  is  to  develop  a  proper 

interpretation  of  the  intrinsic  tbermcdynamic  variables  that  are 

5  14 

conventionally  associate  with  a  given  diatomic  species  ’  and  to 
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delineate  the  role  of  continuum  states  at  high  temperatures. 

We  begin  by  assigning  the  following  partition  function®’^’^®’^^  to 
the  population  of  free  atoms  N^, 


where  is  the  electronic  partition  function  for  the  atom 

lo^Ct)  =  f  cOc 

^  Lzu) 

The  ground  state  is  denoted  by  a=0,  such  that  £^=0.  The  subtleties 

involved  in  the  convergence  of  Eq.  (25),  together  with  the  treatment  of 

17  18 

the  ionization  continua  *  at  high  temperatures  will  be  ignored.  This 
would  merely  distract  us  from  our  primary  goal  of  interpreting  the 

O 

nature  of  A2.  For  the  alkali  atoms,  with  a  doubly  degenerate  5^2 
ground  state,  we  shall  assume  <0^^  =  co^  =  2. 

Note  that  the  entire  volume  V  is  assumed  to  be  accessible  to  the 
center  of  mass  motion  of  the  atom  in  Eq.  (24)  which  we  treat  as  a  point 
particle,  i.e.,  as  an  ideal  ger.  The  total  Hamiltonian  for  the 
interaction  of  tv;o  atoms  A  +  A  is  equivalent  to  the  Hamiltonian  for  the 

-*■  9  1 

diatomic  molecule  A2  for  all  interatomic  distances  R.  ’  Therefore  all 
two-body  imperfections  in  the  vapor  are  contained  in  the  dimer  partition 
function. 


when  the  total  internal  partition  function 


0 


tot 


-  z 


e 


Q,  3  Q: 


(2.7) 


(z$) 
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is  suniiied  over  the  complete  set  of  molecular  electronic  states  i  implied 
by  Eqs.  (l)-(3).  The  subset  of  molecular  states  i(a,a')  which 
dissociate  into  a  particular  pair  of  atoms  A(a)+A(a')  satisfy  the 
constraint, 

.  (30) 


A  pair  of  ground  state  alkali  atoms,  with  a)QWQ=4,  correlates  with  the 
1  +  + 

singly  degenerate  X  and  the  triply  degenerate  b  states  in  Fig.  1. 
These  states  will  be  indexed  as  i=l  and  i=3  respectively  throughout  this 
paper,  and  we  shall  approximate  Q-j-qj  by  Eq.  (29),  consistent  with  our 
assumption  in  Eq.  (26).  Obviously  for  an  accurate  calculation  of  high 
temperature  properties  this  assumption  must  be  relaxed  and  the  complete 
expressions  (25)  and  (28)  must  be  employed. 

For  each  electronic  state  |i>  we  generate  a  set  of 
electronic-rotational  states  li,Jl,mjj^>  with  nuclear  angular  momentum 
quantum  number  t,  and  the  (2t+l )-degenerate  space-fixed  projections  m^. 

A 

These  form  a  complete  set  which  spans  the  polar  coordinates  R  for  the 
orientation  of  the  interatomic  displacement  ^  between  the  interacting 


atoms,  and  we  obtain  an  electronic-rotational  potential  v.  (R)  for  each 

1 

individual  (i,Jl)  state.  For  sufficiently  large  R  the  potentials  will 


o^-oCF?')  -  V;  p  -  ^ 
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where  the  centrifugal  potential  depends  on  t  and  the  reduced  mass  m^/2. 

Often  we  can  approximate  li,t,  mj^>  as  a  simple  product  of  an  electronic 

|i>  and  rotational  Y  (R)  wavefunction 

a.m^ 

U.i.yr.j.'y  «  '^>Y,^Cr) 

with  the  result  that  the  electronic-rotational  potentials  are  additive, 

as  in  Eq.  (6).  In  the  most  general  case,  depending  on  the  particular 
1  9  18  19 

Hund's  coupling  ’  ’  ’  of  electronic-rotational  angular  momentum 

employed  in  v.  -(R),  the  'electronic'  potential  v.{R)  in  Eq.  (6)  may 

1  f  A«  1 

differ  slightly  from  Eq.  (7)  and  exhibit  a  generally  negligible  residual 

dependence  on  A.  However,  the  Hund's  case  (b)  form  in  Eq.  (6)  is 

1  +  3..+ 

perfectly  valid  for  the  E  and  potentials  that  we  shall  ultimately 

y  ^ 

treat,  and  this  form  will  be  employed  throughout  this  study.  Some 

typical  v^  j^(R)  potentials  for  Li2(X'E+)  are  shown  in  Fig.  2. 

The  radial  or  vibrational  motion  of  the  dimer  is  described  by  the 

19  12 

well-behaved  solutions  of  the  Schroedinger  Equation,  ’  * 

['f^^ 

A 

such  that  the  total  electronic-rotational-vibrational  wavefunction 

U./.,rn^.e.y  ■=  (s'/') 

defines  what  we  might  call  the  'internal'  energy  for  the  collision  of 
two  atoms. 

^  +  e  ^35-) 

Below  threshold  for  dissociation  in  the  potential  v.  ,  when  r.<0,  we 


12. 


obtain  a  finite  set  of  true  vibrational'  levels  which  we  index  by 


v=0,l . v„,„(i  ,Ji).These  occur  at  discrete  values  of  E,  such  that 

llidX 


and 


e  = 


hUjcJ) 


(si.b'y 


•<> 


and  the  total  internal  energy  of  the  dimer  is  quantized 


(37) 


Obviously  vibrational  levels  which  lie  above  the  ground  state  limit, 
i.e.,  E.  u>0,  and  asymptotically  correlate  with  excited  state  atoms 

1  j)l>V 
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in  Eq.  (37)  are  subject  to  predissociation  ’  ’  ’  and  the 
definitions  in  Eq.  (36)  are  somev/hat  ambiguous.  However,  the  fact  that 
such  levels  might  not  be  truly  stationary  solutions,  with  discrete 
eigenvalues,  has  a  negligible  effect  on  the  subsequent  evaluation  of 
thermodynamic  properties.  This  phenomenon  is  thoroughly  analyzed  in  the 
next  section  and  should  not  distract  us  from  our  course  at  this  point. 

Above  the  dissociation  limit  the  relative  kinetic  energy  & 
defines  a  continuum  of  wavefunctions  which  exhibit  the  asymptotic 
boundary  condition,^ 


(e.R^ 


A  (3«) 

iv/r  '  2-  '■A*  A  ) 


eO 


where 


(39) 
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The  energy-dependent  elastic  scattering  phase-shift  n.-  .(e)  defines  the 

1 
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elastic-scattering  cross-section  *  for  the  colli  son  of  atoms  A(a)  + 
A(a')  with  relative  kinetic  energy  e,  and  relative  orbital  angular 
momentum,  or  partial  wave  quantum  number  l. 


=  [: 


CJ. 


r«l 


¥ir 
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Well-defined  numerical  procedures  ’  are  available  for  evaluating  the 
phase-shift  from  Eqs.  (33)  and  (38)  as  a  continuous  function  of  e,  and, 
in  particular  we  can  subsequently  evaluate  the  necessary  derivative 
3ni  «/3e  to  arbitrary  accuracy. 

It  should  be  emphasized  that  v;e  have  made  a  severe  approximation  in 

Eqs.  (33)  and  (38)  which  restricts  the  formal  theory.  We  have  neglected 

any  inelastic  scattering  processes  between  the  colliding  atoms  A(o)  + 

A(a')  associated  with  the  electronic  state  i(a,a’)  and  have  assumed  only 

pure  elastic  scattering  is  important.  However,  it  is  our  contention,  based 

on  a  scattering  theory  analysis  of  diatomic  continuum  states  using 

adiabatic  electronic-rotational  (AER)  states,^  that  a  rigorous  treatment 

of  the  inelastic  scattering  will  not  modify  the  thermodynamic 

observables  to  any  significant  degree.  In  any  case,  Eq.  (38)  is 

2  2 

perfectly  valid  for  the  alkali  ground  state  atoms  A(  Sy2)  +  A(  S^^g)  we 

1  +  3  + 

intend  to  study.  There  is  negligible  coupling  between  the  X  Zg  and  b 
continuum  states,  and  we  need  only  consider  the  tv/o  elastic  scattering 
phase  shifts. 
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The  complete  Hamiltonian  for  the  interaction  of  two  atoms  includes 
the  kinetic  energy  operator® associated  with  the  coordinate 
which  locates  the  center-of-mass  position  of  the  collison  complex  with 
total  mass  2m^  within  the  volume  V.  Note  that  the  translational 
partition  function  <ltp(2m^)V  contributed  by  the  center-of-mass 
wavefunction  in  Eq.  (27)  is  approximated  as  a  perfect  gas.  This  is 
valid  if  we  may  assume  that  higher-order  A^+A  and  A2+A2  interactions  are 
negligible  compared  to  A+A  and  do  not  influence  the  activity  of  the 
dimer.  Using  Eq.  (27)  and  (24)  in  (23)  we  obtain  the  equation  of  state 


given  by  Eq.  (15). 


The  internal  partition  function  for  each  i 


th 


electronic  state  in 


Eq.  (28)  is  deconvoluted  into  a  rotational  sum 

^  ^  ^V2.) 

over  the  radial  partition  functions  q^  ,  associated  with  each  specific 
electronic-rotational  potential  v.  «(R)®’®’^,  i.e., 


W-  o,  Ui 
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Equation  (43)  contains  both  bound  q.  ,(B)  and  continuum  q.  (C) 

components  and  measures  the  departure  from  an  ideal  gas.  If  v.=0,  then 

q.  .-0  for  all  A,  and  Q.=0.  This  behavior  is  embodied  in  the  definition 
1  »*>  1 

of  the  phase  shift  n,-  o(e)  in  Eq.  (38).  For  non-interacting  atoms  the 

V.  «  potential  reduces  to  the  simple  centrifugal  potential  shown  in  Fig. 
1 

3  and  the  phase-shift,  which  measures  departure  from  this  ideal 
behavior,  will  vanish  identically.  ’  ’  Thus,  Eq.  (43)  should  be  viewed 
as  the  difference  in  canonical  phase  space  introduced  by  the  interaction 
potential  v^(R).  If  all  interactions  vanish,  then  (p.  =0,  and  we 

retrieve  the  proper  ideal  gas  law  P/kT  =  N/V  for  point  particles  in  Eq. 
(15). 

Each  defined  by  Eqs.  (42)  and  (43)  is  constructed  from  distinct 
bound  (B)  and  continuum  (C)  state  components  as  expressed  in  Eqs. 
(3)-(5),  and  a  comparable  separation  transfers  to  the  total  internal 


partition  function  in  Eqs.  (28)  and  (29), 


Since  the  equilibrium  constant  in  Eqs.  (16)  is  proportional  to 


,2-  « 


(45) 


we  obtain  the  unique  separation  of  and  the  dimer  population 
expressed  in  Eqs.  (20)  and  (21).  The  same  behavior  applies  to  the 


virial  coefficients  in  Eq.  (10),  and  to  the  mean  second-virial 
coefficient  <B>  in  Eq.  (19)  which  is  seen  to  be  a  canonical  average  over 


the  individual  virial  coefficients  associated  with  each  v^(R)  potential. 
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Equation  (47)  applies  to  our  specific  assumptions  for  the  alkali  dimers. 

The  subsequent  breakdown  of  l^gqCC),  Q^-(C),  B^.(C)  and  q^.  (C)  into 

metastable  and  free-continuum  contributions  requires  an  analysis  of  the 

quantum  mechanical  phase-shift  derivative  .(e)/Tr3G. 

1 


III.  CONTINUUM  PHASE  SHIFTS  AND  METASTABLE  STATES 
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The  electronic  potentials  calculated  by  Konowalow^^  for  Li2 
and  Na2  are  shown  in  Fig.  1.  We  shall  use  the  state  of  Li2  to 
demonstrate  the  contribution  of  the  continuum  phase-shifts  and 
metastable  states  to  the  vibrational  partition  function  in  Eq.  (43). 
Specific  electronic-rotational  potentials  v.  (R)  are  plotted  in  Fig.  2 
for  4=0,  100,  130,  and  180  as  the  centrifugal  contribution  to  Eq.  (6) 
progressively  'fills'  the  attractive  v^(R)  potential  until  only  a  pure 
monatonical 1y  repulsive  potential  exists.  Explicit  calculations  will  be 

presented  for  these  representative  angular  momentum  states. 

22  23 

Using  the  Gordon  algorithm  ’  to  solve  Eq.  (33)  we  calculate  that 
a  total  of  v^Q^  =  (Vp,gj^+1)  =  42  bound  vibrational  levels  exist  in  the 
Vi  q(R)  =  v^(R)  potential  which  has  a  dissociation  energy  of  8450  cm"\ 
For  *.=100  the  v^  ^gg  potential  is  only  2560  cm”^  attraction  and  supports 
the  v^p^=ll  true  bound  states  plotted  in  Fig.  4.  However,  behind  the 
rotational  barrier  at  <"0*^  any  t>0  we  may  expect  a  number  of 

quasi-bound  states,  or  metastable  states^’^’^^’^^ to  exist  at  certain 
'eigenvalues'  of  the  incident  kinetic  energy  (See  Fig.  3) 


These  energies  can  be  estimated  by  the  same  procedure  used  to  obtain  the 

true  bound  states  in  Eq.  (35)  by  using  an  approximate  potential  v.  .(R) 

1 

in  place  of  v.  -(R)  in  Eq.  (33),  such  that 
1 


C‘^9') 

. 
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For  £=100  we  obtain  the  ni^jjj.=7  metastable  levels  portrayed  in  Fig.  4. 
We  know  of  course,  that  these  states  are  ultimately  'unstable'  because 

3 

of  tunnelling  through  the  rotational  barrier,  and  we  might  choose  to 
associate  an  energy  width  with  each  level 


Lso) 


which  reflects  the  expected  lifetime  -■  «  „  of  a  molecule  'prepared'  in 

1  5  Jill 

13 

such  a  non-stationary  state-  As  long  as 


JkT  » 


r  • 

‘  A^.JL,yyi 


(S’/) 


we  can  expect  that  the  metastable  nature  of  the  levels  are  of 

insignificant  consequence  and  the  e .  .  „>  0  eigenvalues  can  be  treated 

on  an  equal  footing  with  the  true  bound  states  e.-  <0.  This  can  be 

'  »*>  »v 

demonstrated  by  the  exact  quantum  mechanical  results. 

The  rigorous  solutions  of  Eq.  (33)  for  v^  (R)  yield  a  continuum  of 
wavefunctions  with  the  boundary  conditions  imposed  by  Eq.  (38)  whenever 
e>0.  If  we  were  to  examine  the  amplitude  of  the  wavefunctions  r.  .(e,R) 

If* 

behind  the  barrier  we  v/ould  find  a  negligible  probability  of  observing  a 

pair  of  atoms  with  an  interatomic  displacement  except  in  the 

vicinity  of  es  e.  .  .  This  is  reflected  at  R  -  ~  by  a  rapid  variation 

of  the  phase-shift  -(e)  which  takes  the  following  form,  near  a 

1 
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metastable,  or  resonant  energy,  ’ 


ft-")  » 
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Such  behavior  is  observed  in  the  quantum  mechanical  phase  shift  for 


1=130  in  the  vicinity  of  the  m=8  metastable  level  which  was  calculated 

by  the  procedure  in  Eq.  (49)  to  exist  at  a  "  1947.50  cm  .  The 

24  *  *  -1 

calculated  data  is  plotted  in  Fig.  5.  The  width  r  =  .00255  cm  . 

obtained  from  the  fit  to  Eq.  (52)  corresponds  to  a  lifetime  of  10"^  sec, 

and  the  position,  «  _  =  1947.775  cm"\  is  only  slightly  shifted  from 

the  approximate  result.  Note  that  the  phase  increases  by  modular-ir  in 

passing  from  e<e,.  .  _  to  e>e,-  .  Figure  6  shows  similar  behavior  for 

the  top-most  resonance,  m=9,  which  lies  about  15  cm”^  below  the 

rotational  barrier  at  =  2056  cm"\  In  this  case  the  width  obtained 

-1  -1 2 
from  Eq.  (52)  is  r23.2  cm“  ,  with  a  resultant  lifetime  of  10  sec. 

The  derivative  3n,-  a/ttSe  simply  yields  the  normalized  Lorentzian 

I  fX 
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line  shape  *  ’  associated  with  a  width  r.  .  ,  i.e.. 


TT  'a  e 


/  ztr 


re-e. 


1st) 


and  the  contribution  to  the  continuum  partition  function  (31)  in  the 
vicinity  of  any  resonance  is  bounded  by  unity,  i.e.. 


I 


c/e  e  -3^/./. 


"A 


rr'6e. 


Note  that  similar  considerations  will  apply  to  any  metastable  or 
predissociating  state  contributing  to  the  dimer  partition  functions,  and 
not  merely  to  rotational  predissociation  through  centrifugal  barriers; 
the  only  restriction  to  consider  is  the  criterion  in  Eq.  (51).  It 
should  also  be  emphasized  that  the  role  of  niatastable  molecules  in 


equilibrium  properties  is  independent  of  the  lifetime  of  the  state,  and 

3  4 

is  not  contingent  on  any  dynamic  assumptions,  ’  such  as  requiring  the 

collision  frequency  between  free  atoms  to  exceed  the  predissociation 

rate.  This  can  be  appreciated  by  examining  the  nature  of  the 

phase-shift  more  closely. 

21 

Scattering  theory  predicts  that  two  colliding  atoms  A{a^)+A(a) 
progressing  on  the  molecular  potential  v^{R)  with  orbital  angular 
momentum  Jl,  and  relative  kinetic  energy  e,  will  experience  a  'delay 
time'  or  collision  lifetime  (At)^  ^  ^  relative  to  two  non-interacting 
atoms  with  the  same  dynamic  variables,  such  that 


Thus  the  continuum  partition  function  q.  »(C)  in  Eq.  (43)  may  be 

1 
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expressed  as  follows, 

^£/JkT  ^  ^ 

CL:  g  (C^  =  J  ^ 

o  h 


A  positive  delay  time  means  a  positive  contribution  to  the  partition 
function  and  a  higher  density  of  states  compared  to  the  free  atom 
interactions.  This  is  particularly  evident  in  the  vicinity  of  a 


resonance  where  Eq.  (53)  predicts  a  peaked,  positive  delay  time 


(At) 


which  increases  with  the  lifetime  of  the  metastable  state.  In 


(^7) 


contrast,  a  negative  delay  time,  as  occurs  for  a  repulsive  hard-core 

potential,  imposes  a  negative  contribution  on  q.  .(C),  which  is,  in 

1 

fact,  equivalent  to  the  net  loss  of  pfase  space  associated  with  the 


excluded  volume  defined  by  the  interaction  potential. 


For  real  potentials,  such  as  in  Fig.  2,  there  v/ill  be  both  negative 
and  positive  contributions  due  to  the  repulsive  and  attractive  regions 
of  interaction.  However,  we  may  use  the  analysis  in  Eqs.  (53)  and  (54) 
to  isolate  the  dominant  positive  contribution  of  the  resonance  states  to 
the  continuum  partition  function,  i.e.. 


m 


O0 

de 


^  -fr**- 
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The  'free-continuum'  component  q.  ^^(F)  embodies  all  the  remaining 
non-resonant  energy  dependent  contributions  of  ri<  A(e)  between  the 

1  >jC 

resonances  for  in  Fig.  4,  as  well  as  the  entire  j^(e) 

3  * 

dependence  above  the  rotational  barrier.  Unless  one  is  devising 
experiments  in  which  the  dimer  metastable  states  are  directly 
'observable'  the  only  pertinent  quantity  in  the  analysis  of  the 
thermodynamic  properties  is  the  total  continuum  partition  function  Q^(C) 
and  this  deconvolution  into  components  is  only  for  our  convenience. 

If  the  interaction  v^(R)  between  two  atoms  vanishes  only  continuum 
solutions  to  Eq.  (33)  exist  and  the  asymptotic  form  of  the  radial 
function  in  Eq.  (38)  is  equivalent  to  the  spherical  Bessel  function 
1.8,9,12,21  j  (kf.)  for  the  partial  wave  t,  where  the  phase  shift  is 
identically  zero.  Thus  the  phase-shift  for  an  arbitrary  potential  is  a 
measure  of  the  deviation  from  the  ideal  gas,  or  point  particle, 
interaction  which  leads  to  the  definition  of  the  delay  time  in  Eq.  (55). 

The  boundary  condition  in  Fq.  {'.6)  is  indeterminate  by  modular  tt. 


and  we  must  convert  the  principal  value  of  n.-  obtained  from  the 
numerical  solution  of  Eq.  (33)  into  a  continuous,  analytic  function  of 
e>0  which  can  be  differentiated  and  used  in  Eq.  (56).  Levinson's 
theorem^ demonstrates  that  the  change  in  the  phase  shift 
between  e=0  and  e=<»  is  simply  related  to  the  total  number  of  bound 
states  Vtott"*  »*')  exist  in  the  v^  ^^^(R)  potential, 

and,  since  n-  ,(“)=0  for  molecular  potentials  v.  which  are  typically 

1  )  X  1 

finite  at  R=0,  we  require  a  threshold  value  of  ^(0)  =  ^^^^^(1,*-). 

The  actual  threshold  behavior  for  a  potential  which  approachs  v^-  -  CR“^ 
at  large  R  is  known  to 

•7^*  -fori  ^0,1 

£-=>  O 

~  *0-(k‘*)  4  ov  ^  2  (i>oV 

and  we  could  analytically  continue  the  numerically  derived  phase  shifts 
to  obtain  a  unique  value  of  n(e)  for  all  e. 

Fortunately  we  have  found  an  excellent  approximation,  based  on  the 

c  no 

WKB  approximation,  ’  which  supplies  a  simple,  analytic,  single-valued 

estimate  of  n-  .(c)  at  any  energy. 

1  *x 


The  symbol 


Jo 


s  meant  it  indicate  the  exclusion  of  all  non-classical  . 


regions  of  phase  space  from  the  integrals,  which  occurs  v/henever  V 
becomes  imaginary.  When  used  in  Eq.  (56)  and  summed  over  I  in  Eq. 

(42),  this  approximation  is  equivalent  to  the  classical  contribution  of 
B^(C)  to  Eqs.  (12)  and  (13).  In  the  following  plots  we  will  compare 
both  n,-  «  and  and  find  extraordinary  agreement  except  right  at 
threshold,  and  in  the  vicinity  of  resonances.  As  already  seen  in  Figs. 

5  and  6,  Eq.  (61)  transforms  the  step-like  increases  in  j^(e)  at 
e  ~  e,-  0  m  3  smoothly  varying  function.  This  is  equivalent  to 
reversing  the  procedure  in  Eq.  (54)  and  substituting  a  continuous 
function  to  represent  the  resonant  contribution  to  the  partition 
functions.  The  final  results  will  be  equivalent  whenever  kT  > (e^  Jl  m+1 
-  quantitively  justify  the  use  of  the  classical 

theory. 

The  phase-shift  for  1=0  is  plotted  in  Fig.  7.  It  falls  fairly 
rapidly  from  its  threshold  value  of  42Tr,  introducing  a  large  negative 
contribution  to  Eq.  (56).  At  extremely  large  energies  the  phase 
actually  becomes  negative  as  the  repulse  core  of  the  v^  ^  potential  in 
Fig.  2  becomes  dominant,  and  then  slowly  progresses  to  zero  as  A 

simple  physical  interpretation  of  the  negative  contribution  near 
threshold  is  that  the  rapid  acceleration  of  the  atoms  as  they  pass  over 
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the  deeply  attractive  well  causes  the  atoms  to  spend  less  time  in  this 

region  of  phase  space,  i.e.,  (At).  .  <0,  compared  to  free  atoms.  It  is 

•  %x>  »e 

Interesting  to  note  that,  at  infinite  temperature,  there  is  an  exact 
cancellation  of  the  positive  bound  state  contribution 

^Ttot 

ot> 

and  the  negative  continuum  state  contribution  predicted  by  Levinson's 
theorem, 

r-"  ■  ‘^tot 

such  that  q.  in  Eq.  (43)  is  identically  zero.  In  any  case  it  should 

•  jJC 

be  apparent  that  neglect  of  the  continuum  phase  space  in  estimating  q. 

1  ,0 

at  high  temperature  will  overestimate  the  magnitude  of  the  partition 

function.  This  is  demonstrated  in  Fig.  8  where  the  various 

contributions  of  q.  .(B),  q.  .(H)  and  q,.  .(F)  for  Li-,(x\'^)  are  plotted 
1  I  ^  y 

as  a  function  of  t  for  T  =  6000®K.  A  similar  plot  for  ^^(^E^)  is  shown 

z'  u 

in  Fig.  9.  These  quantities  v/ere  actually  calculated  using  the  WKB 
theory  but  the  results  are  in  quantitative  agreement  with  the  quantum 
values. 

It  should  be  pointed  out  that,  for  homonuclear  molecules  with 
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nuclear  spin  1=0,  only  even  values  of  H  satisfy  the  antisymmetrization 
of  the  total  dimer  v/avefunction,  and  alternate  values  of  t  are  missing 
from  the  summation  in  Eq.  (42).  This  effect  is  equivalent  to  the 
ubiquitous  'symmetry  factor'  a=2  which  is  used  to  reduce  the  magnitude 
of  the  classical  partition  functions  for  homonuclear  species. 

The  fact  that  Li  has  a  nuclear  spin  1=1  gives  rise  to  the  same 
complications  that  occur  in  the  D2(X^Zg)  state^^’^^  with  an  alternating 
weighting  of  2/3:l/3  assigned  to  Jl=even  and  ji=odd  eigenvalues 
respectively.  However,  the  evaluation  of  statistical  properties  at  high 
temperatures  is  quite  insensitive  to  this  subtlety  and  we  shall  simply 
sum  over  even  values  of  i  with  a  weighting  of  unity. 

The  phase  shift  for  £=100  is  shown  in  Fig.  10  and  should  be 
examined  in  comparison  to  Fig.  4.  The  threshold  value  Hit  reflects  the 
presence  of  v^q^=11  bound  states  in  the  v^  loo^^^  potential.  There  is  a 
rapid  increase  in  the  phase  to  a  value 

rotational  barrier  at  e^slOOO  cm’*^  concomitant  with  the  7  metastable 
states  which  are  predicted  to  exist  behind  the  barrier.  The  WKB  phase 
shift  predicts  a  smooth  monotomic  increase  from  threshold  to  e  =  1000 
cm”\  The  actual  quantum  mechanical  phase  shift  behaves  as  in  Fig.  5 
and  remains  relatively  flat,  with  an  abrupt  series  of  step-like 
increases  by  modular-tr  at  each  metastable  resonance  position.  Above  e  = 
1000  cm“^  both  phases  coincide  and  contribute  first  positive,  < 
0,  and  their  negative, >  0,  components  to  the  integral  in  Eq.  (56). 
The  resultant  q^  loo^*'^  ^i  100^^^  contributions  at  T=6000°  can  be 
seen  in  Fig.  8. 

The  phase  for  £=130  in  Fig.  11  must  obviously  start  at  «(0)=0 
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since  the  potential  in  Fig.  3  is  seen  to  support  no  bound  states  below 
threshold.  There  is  a  negligible  increase  in  phase  until  the  minimum  in 
the  12Q  potential  becomes  accessible  at  e=s:706  cm”^  and  then  there  is 
a  rapid  rise  to  Qtt  as  the  nine  metastable  states  contribute  to  the 
scattering.  Again  we  have  plotted  the  smoothly  varying  WKB  phase  which 
smooths  the  step-like  quantum  phase.  We  shall  argue  that,  for  kT  large 
compared  to  the  metastable  spacings,  the  monotonic  function,  which  leads 
to  exact  agreement  with  the  classical  partition  function,  is  an  adequate 
approximation. 

Finally,  we  present  the  £=180  phase  in  Fig.  12.  As  seen  in  Fig.  2 
this  potential  supports  neither  bound,  nor  metastable  states,  and  we 
obtain  a  smoothly  varying  phase-shift  which  predominately  makes  positive 
contributions  to  the  partition  function  due  to  the  deeply  attractive 
well  in  the  v.(R)  electronic  potential. 

The  calculated  electronic  partition  functions  for  the  singlet,  Q^, 
and  triplet,  Q^,  states  of  Li2  and  Na2  are  presented  in  Tables  I  and  II. 
The  tabulated  ratios  Q^(M)/Q^{B)  and  Q^(F)/Q^(B)  give  a  measure  of  the 
importance  of  the  various  components  of  the  continuum  states  as  the 
temperature  is  increased.  Obviously  these  dissociative  states  are  most 
significant  in  the  shallow  b  potential,  but  the  neglect  of  Q^(M),  or 
Q^(F)  can  not  be  justified  even  for  the  strong,  chemically  bonded  X^E^ 


state. 


IV.  DISCUSSION  OF  RESULTS 


Given  the  partition  functions  for  Lig  and  Nag  we  can  derive  the 
usual  thermodynamic  variables  for  these  molecules.  We  have  compared  our 

5 

calculations  to  the  JANAF  tables  and  find  substantial  deviations  at 
temperatures  above  2000°K.  We  concluded  that  any  tabulations  of 
thermodynamic  properties,  such  as  the  JANAF  tables,  which  rely  on  a 
Dunham  expansion  of  the  molecular  term  values  must  ultimately  fail  at 
high  temperatures.  In  particular  the  anharmonicity  substitutes  a 
spuriously  high  density  of  discrete  states  in  place  of  the  true 
continuum  states.  These  meaningless  states  tend  to  over  estimate  the 
molar  entropy  and  enthalpy.  This  fault  is  especially  evident  in  the 
tabulated  molar  heat  capacity  Cp  which  often  increases  monotonically  to 
unrealisticly  high  values.  More  sophisticated  calculations^®  which 
attempt  to  evaluate  the  finite  sums  over  bound,  Eq.  (4),  and  metastable, 
Eq.  (9),  states  yield  substantially  more  reliable  results,  but  this 
approach  still  suffers  from  the  neglect  of  the  free-continuum 
contribution  Q^.(F)  in  Eq.  (8).  This  can  be  seen  in  Tables  I  and  II. 

In  this  discussion  we  shall  consider  the  one  observable  parameter 
that  is  derivable  from  the  thermodynamic  analysis;  the  equilibrium 
constant  or  equivalently  the  second  virial  coefficient  <B>  = 

We  can  not  measure  the  molar  Gibb's  free  energy,  or  enthalpy,  or  entropy 
of  an  isolated  alkali  dimer.  Hov/ever  we  can  confine  a  prescribed  number 
of  atoms  N  in  a  volume  V  at  temperature  T  and  measure  the  total  pressure 
P  of  the  resultant  vapor.  Thus  we  car,  obtain  the  total  population  of 
dimers  from  the  equation  of  state  (15)  and  the  constraint  in  Eq.  (14) 


^AA  =  '^A«  ( =  "PaaV/ji^t  =(N-TV/f?T).  ((,4) 


Given  our  definition  of  the  dimer  embodied  in  our  construction  of  the 
partition  functions  (24)  and  (27)  we  may  equate  in  Eq.  (16)  to 
directly  measurable  quantities. 


Ke,  =  hr  -  AT  Ta« 


(6S) 


The  various  components  of  in  Eq.  (20)  manifest  themselves  directly 
in  the  predicted  dimer  population  (21)  and  make  explicit  contributions 
to  the  observable  equation  of  state  (22). 

The  solid  curves  in  Fig.  13  present  for  Li2  and  Na2  when  the 
entire  contribution  of  singlet  and  triplet  ground  states  are  included  in 
Eq.  (47).  These  are  calculated  using  the  classical  expression  in  Eq. 
(12).  This  is  the  equilibrium  constant  v/e  would  expect  to  extract  from 
Eq.  (65)  assuming  that  no  excited  electronic  states  of  the  atoms  and 
molecules  influence  our  observations  of  in  Eq.  (64).  Actually 
electronic  excitations  seem  to  influence  the  thermodynamics  above  about 
4000°K^^’^^  and  our  results  at  higher  temperatures  must  be  viewed  as 
only  a  partial  analysis  of  the  real  alkali  vapor.  The  calculations  in 
the  range  1000-3000°K  should  be  quite  reliable  and  are  summarized  in 


Table  III. 

The  dashed  curve  in  Fig.  13  is  the  resultant  K  obtained  if 

3 

we  exclude  the  contribution  in  Eq.  (47).  In  this  case  Eq.  (16) 

only  yields  Ny^yy  =  Ny^y^  Oz*)  and,  analogous  to  Eq.  (22),  imperfections 

^  3+2 

due  to  the  missing  triplet  contribution  fly^y^  (  l^)  =  -3B2Ny^/(4V)  must 
yield  non-ideal  behavior  in  the  monomer  component  of  the  equation  of 
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state.  The  effect  is  particularly  pronounced  in  the  Na2  vapor,  where 
the  large  negative  contribution  due  to  the  free-continuum  component 
Q3(F)  in  Table  II  causes  a  48%  reduction  in  at  10,000°K. 

Fortuitously  this  just  about  cancels  the  positive  Q^(F)  contribution  due 
to  the  attractive  singlet  state  and  it  might  seem  reasonable  to  use  the 
usual  approximation  ’  *  Q-|(M)  and  neglect  all 

free-continuum  contributions  entirely.  However  in  Table  I  we  see  that 
Q3(F)  is  much  reduced  in  Li^  and  there  is  never  any  assurance  that  a 
truncated  summation  over  bound,  or  bound  and  metastable  states,  will 
lead  to  an  accurate  The  only  rigorous  procedure  is  to  completely 

evaluate  both  and  Q^. 

Direct  comparisons  between  this  calculation  and  the  Kp  tabulated  in 
14 

the  JANAF  Tables,  which  employ  Dunham  type  expansions  for  the 
molecular  eigenvalues,  are  complicated  by  many  factors,  and  especially 

by  the  substantial  discrepancy  between  the  H°290  predicted  by  the 
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Konowalow  et  al  potentials,  and  the  old  values  of  H  293  employed  in 

the  Tables.  Since  this  quantity  is  systematically  subtracted  from  the 

tabulated  molar  thermodynamic  variables  for  the  dimer,  comparisons  with 

the  JANAF  Tables  are  best  considered  in  reference  5. 

Obviously  the  thermodynamics  of  Li2  is  dominated  by  the  deeply 
attractive  state  in  Fig.  1,  with  D^  =  8450  cm"\  If  we  ignore  this 
state,  and  only  consider  the  equilibrium  properties  associated  with  the 
weakly  attractive  potential  (D^  =  292  cm  ),  we  can  simulate  the 
behavior  of  a  loosely  bound  dimer.  This  is  done  in  Fig..  14,  where  the 
solid  curve  gives  the  equilibrium  constant  obtained  by  using  <B>  =  in 
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Eq.  (47).  The  dashed  curve  results  front  excluding  the  continuum,  and 
using  <B>  =  82(6).  The  qualitative  features  of  these  calculations  can 
be  expected  to  mimic  the  equilibrium  properties  of,  say,  the  group  IIA 
dimers,  Be2,  Mg2,  Ca2.  and  the  group  I  IB  dimers,  Zn2,  Cd2,  Hg2,  of  the 
periodic  table  all  of  which  have  shallow  van-der-Waal  ground 
states. 

Because  of  the  dominant  contribution  of  the  free-continuum 
component  B2(F),  the  excluded  volume  introduced  by  the  repulsive  portion 
of  the  potential  causes  Bo  to  become  positive  at  5400°K.  This 
result  could  be  anticipated  from  Fig.  8,  where  the  individual  rotational 
contributions  to  are  negative  for  all  a.  For  such  a  system  the  molar 
thermodynamic  properties,  such  as  6°  “  lnQ3»  ill  defined.  We  could 
content  ourselves  by  defining  the  dimer  in  terms  of  the  bound-state 
portion  of  the  partition  function  Q^CB)  in  Tables  I  and  II,  or  possibly 
by  the  combination,  Q3(B)  +  Q^CM),  and  tabulate  the  resultant 
thermodynamic  variables  accordingly.  However,  without  an  understanding 
that  the  'monomer'  imperfections  attributable  to  B^CF)  make  a  dominant 
contribution  to  Eq.  (22b)  we  could  not  hope  to  use  the  thermodynamic 
tables  to  predict  the  observable  properties  of  the  vapor. 

After  this  elaborate  analysis  of  diatomic  thermodynamics,  let  us 
ask  what  constraints  we  should  place  on  the  use  of  current  thermodynamic 
tabulations,  and  what  aspects  of  this  polemic  can  be  ignored? 

Obviously  at  low  to  intermediate  temperatures,  when  Q^q^(B)» 
Qtot^C)’  we  can  ignore  the  entire  continuum  and  metastable  problem.  At 
all  pressures  the  conventional  spectroscopic  analysis,  using  the  Dunham 
expansion,  should  be  adequate. 

-  31  - 


At  higher  temperatures  the  importance  of  the  continuum  states 

depends  on  the  ratio  elevated  pressures,  with  large 

departures  from  the  ideal  gas  law  PV  =  NkT,  the  deviations  from  the 

JANAF  predictions  can  be  substantial.  Particularly  if  our  goal  is  to 

obtain  a  quantitative  description  of  the  equation  of  state  our  only 

course  is  to  follow  the  prescriptions  in  this  paper  and  evaluate 

accurate  partition  functions  using  Eq.  (12).  Of  course  in  this  case 

the  analyses  must  be  extended  to  include  assessments  of  higher  order 
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interactions  in  the  vapor.  For  instance,  it  is  known  that  the  alkali 
trimers  are  significantly  bound.  However,  this  only  emphasizes  the  need 
for  accurate  evaluation  of  the  dimer  equilibrium,  especially  if  we  hope 
to  use  apparent  departures  from  the  expansion  in  Eq.  (15)  to  estimate 
these  higher-order  effects. 

At  high  temperatures  and  modest  pressures,  when  2N^^/N^  is 
appreciably  less  than  one,  we  would  still  require  accurate  partition 
functions  to  obtain  a  reliable  equation  of  state.  However,  if  the 
thermodynamics  is  to  be  used  merely  to  predict  concentrations  of 
specific  internal  quantum  states  A2(i,t,v)  with  electronic-rotational- 
vibrational  energy  E.  i.e., 

*  *  C"  '  /CtoT 

small  deviations  of  =  N  [1-2N^/N]  from  N  may  often  be  ignored, 
or  adequately  corrected  using  existing  thermodynamic  data.  However,  we 
now  require  accurate  evaluation  of  the  molecular  eigenvalues  E.  .  ,  and 

1  >Xr  J  V 

at  high  temperatures  we  can  expect  an  appreciable  percentage  of  dimer 
molecules  to  exist  in  the  continuum  states  associated  with  the  partition 
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functions  ^(C)  and  Q^(C)  in  Eqs.  (43)  and  (5).  In  the  case  of 

3 

van-der-Waals  type  molecules,  such  as  the  state  of  Li2,  this 
temperature  is  not  remarkably  high,  and,  as  seen  in  Table  I,  even 
for  the  deeply  attractive  ground  state  of  Li2  approximately  15%  of 
the  molecules  are  metastable,  and  12%  exist  in  the  free-continuum  at 
6000°K. 
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V.  SUMMARY  AMD  RECOMMENDATIONS 


At  high  temperatures  the  metastable  state  partition  function, 
Qtot(M)»  the  free-continuum  partition  function,  Q^q^(F')» 
important  components  in  the  total  internal  partition  function,  = 
Qtot(B)  ^tot^*^'^  ^  ^tot^*^^’  ®  diatomic  molecule.  We  could  choose 

to  represent  the  dimer  simply  in  terms  of  the  finite  set  of  true 
molecular  bound  states  which  define  but  then  we  must  be 

cognizant  of  the  fact  that  all  continuum  contributions  to  the  molecular 
phase  space,  merely  been  transferred 

to  a  new  definition  of  the  monomer  which  can  cause  marked  imperfections 
in  the  activity  of  what  we  now  choose  to  call  the  atom  population. 

Although  the  metastable  component  "is  always  positive,  the 

remaining  contributions  of  the  continuum,  Q^q^(F)»  can  be  either 
positive  or  negative  with  a  magnitude  that  can  match  or  exceed 
Thus  there  is  no  assurance  that  a  dimer  defined  by  the  combination  of 
bound  and  metastable  states,  Q^p^(B)  +  will  yield  any  superior 

results. 

The  safest  course  is  to  use  the  total  partition  function  to 
define  the  thermodynamic  properties  of  the  dimer  since  there  is  never 
any  doubt  that  the  monomer  is  properly  described  as  a  perfect  gas.  In 
any  case  it  is  imperative  to  use  a  mutually  consistent  scheme  when 
defining  the  dimer  partition  function  for  diverse  species.  Otherwise 
substantial  errors  can  be  introduced  into  the  calculation  of  equilibrium 
distributions  in  mixed  vapors. 
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We  offer  two  recommendations  for  the  calculation  of  thermodynamic 
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properties  at  high  temperatures: 

1)  All  high  temperature  calculations  should  begin  with  the  interatomic 

potentials  v^(R)  for  each  electronic  state  included  in  the  partition 

function.  The  spectroscopic  parameters  obtained  from  an  expansion  of 
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observed  term  values  can  be  used  as  input  to,  say,  an  RKR  potential 
for  the  state.  However,  it  is  imperative  to  use  potentials,  such  as  in 
Fig.  2,  which  properly  represent  the  dissociation  process  and  therefore 
introduce  the  correct  density  of  states  for  the  molecular  continue. 

2)  We  have  shown^  that  the  classical  approximation  to  the  partition 
function  is  extremely  accurate  at  high  temperatures.  The  simplest 
course  is  to  use  the  electronic  potentials  v^(R)  in  the  classical 
expressions  and  numerically  evaluate  the  complete  partition  function 
in  Eq.  (12).  This  result  contains  all  contributions  to  the  molecular 
phase  space,  including  the  metastable  and  free-continuum  components. 

The  monomer  is  then  rigorously  represented  as  an  ideal  gas,  and  all 
imperfections  are  embodied  in  our  definition  of  the  dimer.  For  shallow 
van-der-Waal  potentials,  where  may  become  negative,  there  will  be 
some  subtleties  in  this  procedure.  However,  the  analysis  we  have 

5 

presented  in  terms  of  virial  coefficients  for  the  monomer  can 
rigorously  handle  this  occasional  difficulty. 


ACKNOWLEDGEMENT 

This  work  is  supported  in  part  by  the  Air  Force  Office  at 


Scientific  Research,  Contract  Nc.  ArGSR-ISSA-80-0012. 


References 


1.  F.  H.  Mies,  Molec.  Phys.  41,  973  (1980) 

2.  T,  L.  Hill,  J.  Chem.  Phys.  n,  617  (1955) 

3.  D.  E.  Stogryn  and  J.  0.  Hirschfelder,  0.  Chen.  Phys.  31,  1531 

(1959) 

S.  K.  Kim  and  J.  Ross,  J.  Chem.  Phys.  42,  263  (1965) 

5.  F.  H.  Mies  and  P.  S.  Julienne,  submitted  J.  Res.  Nat.  Bur.  Stand. 

6.  0.  K.  Rice,  Statistical  Mechanics  Thermodynamics  and  Kinetics 

(W.  H.  Freeman,  San  Francisco,  1967)  chapters  8,  11,  16. 

7.  J.  0.  Hirschfelder,  C.  F.  Curtiss  and  R.  B.  Bird,  Molecular  Theory 
of  Gases  and  Liquids  (Wiley,  Mev/  York,  1954)  chapter  6 

8.  N.  F.  Mott  and  H.  S.  W.  Massey,  The  Theory  of  Atomic  Collisions 
(3rd  Ed.,  Oxford,  London,  1965) 

9.  F.  H.  Mies,  Molec.  Phys.  H,  953  (1980) 

10.  R.  C.  Feber  and  C.  C.  Herrick,  U.S.  AEC  LA  3597  (1966) 


11.  L.  V.  Gurvich,  I.  V.  Veits,  V.  A.  Medvedev,  G.  A.  Khachkuruzov, 

V.  S.  Yungman,  G.  A.  Bergman  et  al,  "Termodinamicheskie  Svoistva 
Individual'nykh  Veshchestv",  (Thermodynamic  Properties  of 
Individual  Substances),  V.  P.  Glushko,  gen.  ed.,  Vol.  1,  part  1 
(1973)  Izdatel'stvo  "Nauka"  Moscow 

12.  F.  H.  Mies  arid  P.  S.  Julienne,  in  preparation 

13.  U.  Fano,  Phys.  Rev.  1^,  1866  (1961) 

14.  D.  R.  Stull  and  H.  Prophet,  JAflAF  Thermodynamic  Tables,  (2nd  Ed., 
Nat.  Stand.  Ref.  Data  Ser. ,  Nat.  Bur.  Standards  37,  1?71 ) 

15.  D.  0.  Konowalow  and  M.  L.  Olson,  J.  Chem.  Phys.  450  (1979); 

D.  D.  Konowalow,  M.  E.  Rosenkrantz  and  M  L.  Olson,  J.  Chem.  Phvs. 
72,  2612  (1980) 

16.  G.  S.  Rushbrooke,  Introduction  to  Statistical  Mechanics  (Oxford, 
London,  1949) 

17.  K.  S.  Pitzer,  J.  Chem.  Phys.  70,  393  (1979) 


■  36  " 


18.  L.  D.  Volyak,  Thermodynamic  Properties  of  Gases  and  Liquids, 

V.  A.  Rabinovich,  Ed.,  (Committee  of  Standards,  USSR,  Translated  by 
Israel  Program  for  Scientific  Translations,  1970,  TT69-55091)  p.  43 

19.  b.  Herzberg,  Spectra  of  Diatomic  Molecules,  (2nd  Ed.,  Van  Nostrand, 
Princeton,  1950] 

20.  P.  S.  Julienne,  J.  Molec.  Spectros.  270  (1975) 

21.  T.-Y.  Wu  and  T.  Ohmura,  Quantum  Theory  of  Scattering 
(Prentice-Hall,  Englewood  Cliffs,  N.J.  1962) 

22.  R.  Gordon,  J.  Chem.  Phys.  51,  14  (1969);  Methods  Comput.  Phys.  !£, 
81  (1971) 

23.  R.  I.  LeRoy  and  R.  B.  Bernstein,  J.  Chem.  Phys.  4312  (1968); 

R.  B.  Bernstein,  J.  Chem.  Phys.  33,  795  (1960) 

24.  We  chose  a  principal  value  for  tan”^X  such  that  tan~^(-“)=  0  and 
tan"  (+«>)  =-17 

25.  This  quantum  expression,  which  replaces  each  microcanonical  volume 
of  phase  space  with  an  average  time-delay  for  each  energy  state, 
bears  some  interesting  comparisons  to  the  ergodic  and  semi-orgodic 
problem;  e.g.,  R.  C.  Tolman,  The  Principles  of  Statical  Mechanics 
(Oxford,  London,  1938);  T.  L.  Hill,  Statistical  Mechanics 
(McGraw-Hill,  New  York,  1956) 

26.  J.  R.  Taylor,  Scattering  Theory:  The  Quantum  Theory  of 
Nonrelativistic  Collisions  (Wiley,  fiew  York,  1972) 


27.  M.  L.  Goldberg  and  K.  M.  Watson  Collision  Theory  (Wiley,  New  York, 
1964) 

28.  L.  D.  Landau  and  E.  H.  Lifshitz,  Quantum  Mechanics: 
Non-Relativistic  Theory  (Pergamon  Press,  London,  1958) 

29.  F.  W.  Reiter,  Z.  Naturforsch.  28a ,  1676  (1973) 

30.  R.  I.  Artyin,  High  Temp.  6,  966  (1968);  translated  from  Teplofizika 
Vysokikh  Temp.  6,  11010  X1968) 

31.  C.  H.  Wu,  J.  Chem.  Phys.  3181  (1976);  R.  L.  Martin  and 
E.  R.  Davidson,  Molec.  Phys.  1713  (1978) 

32.  R.  N.  Zare,  A.  L.  Schmeltokopf ,  H.  J.  Harrof,  and  D.  L.  Albritton, 
J.  Molec.  Spectros.  37  (1973) 

34 

33.  These  rcconrapndations  were  originally  proposed  by  C.  W.  Beckett 
and  L.  Haaf^'’,  but  ha  not  been  properly  exploited. 


-  37  - 


34.  C.  W.  Beckett,  and  L.  Haar,  Proc.  of  the  Joint  Conference  on 
Thermodynamic  and  Transport  Properties  of  Fluids,  Inst,  of 
Mechanical  Engineers,  London,  p.  27  (1957). 

35.  L.  Haar,  "Kinetics,  Equilibria  and  Performance  of  High  Temperature 
Systems,"  Proc.  of  the  First  Conference  of  the  Western  States 
Section  of  the  Combustion  Institute,  Butterworth  and  Co.,  Limited, 
London,  1959. 


-  38  - 


(b)  units  cal.  mole”  deg."  .  Evaluated  using  H2gg  =  -21.480  kcal.  and  assuming  =  Q-]  +  SQg 


TABLE  III.  Equilibrium  Constant.  Enthalpy  and  Gibbs  Energy  of  Formation  for  LI.  and  Na. 
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Figure  Captions 


Fig.  1.  Potentials  for  various  Li^  and  Mag  electronic  states 
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calculated  by  Konowalow  et  al. 

Fig.  2.  Representative  electronic-rotational  potentials,  v.  «(R),  for 

1 

the  state  of  Li«. 

g  2 

Fig.  3.  Schematic  comparison  of  a  typical  electronic-rotational 

potential  v^  ^(R)  and  a  pure  centrifugal  potential  Yjj_(R)  = 

hH(£+l)/m^R^  associated  with  noninteracting,  i.e.,  ideal, 

atoms.  The  phase-shift  n-  n(e)  for  a  pair  of  incident  atoms 

1  ,* 

with  relative  kinetic  energy  e  measures  the  departure  from 

ideal  gas  scattering.  A  positive  n..  a  is  associated  with  an 

1 

increase  in  the  number  of  nodes  in  the  scattering 
wavefunction  w'l^never  v^  is  more  attractive  than  At 
energies  above  the  crossing  of  these  two  curves  we  expect  a 
negative  phase-shift  due  to  the  exclusion  of  the  interacting 
atoms  from  the  phase  space  available  to  the  ideal  atoms. 

Fig.  4.  The  calculated  bound,  ^  <  0.  and  metastable. 

Cl  ^QQ  >  0,  states  are  located  in  a  plot  of  v^ 
the  £=100  state  of  Li2(X^z:g).  The  vibrational  levels  v=0,10 
define  the  true  bound  state  portion  radial 

partition  function  "  '^1,100^^^  ‘’l,100^^^’ 

continuum  contribution  q^^go^^^  "  ‘^1,100^'^^  ^1,00^*"^ 
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divided  into  two  components.  The  finite  set  of  virtual  levels 
in=l ,7  predicted  to  exist  behind  the  rotational  barrier 
R  <  ,  and  with  approximate  eigenvalues  0<e^  ^qq 

,  define  the  metastable  contribution 
remainder  of  the  phase-space  associated  with  regions  to  the 
right,  and  above  the  rotational  barrier  is  assigned  to  the 
free-continuum  component 


Fig.  5.  Elastic  scattering  phase-shift  vicinity  of 

the  predicted  resonance  at  130  m  ~  1949.50  cm"^  for  the 

£=130  rotational  state  of  Li2(X^Zg).  This  is  the  8th 

resonance,  m=8,  which  lies  about  53  cm”  below  the  rotational 

barrier  for  v^  130^^^  figure  presents  the 

exact  quantum  mechanical  predictions  for  analytic 

fit  to  n  =  (background)  -  tan”^  P  jL  r  /(e-£„)^  perfectly 

reproduces  the  numerical  results.  The  extracted  width 

^1  130  8  ”  0*0051  cm”^  predicts  a  predissociation  lifetime 
=  10-9 

Ti  120  g  sec.  for  this  particular  metastable  state. 


Fig.  6.  Elastic  scattering  phase  shift  ni  i3q(^)  in  the  vicinity  of 

the  the  top-most  resonance,  m=9,  for  £=130  rotational  state  of 

Li2(X  Eg)*  Such  resonances  which  lie  close  to  the  rotational 

barrier  (2057  cm”^)  are  often  broad,  e.g.,  r«3  cm”^  and 
12 

T  «  10  sec.,  and  are  not  well  fitted  to  the  simple 
expression  used  in  Fig.  5.  However  the  width  is  still  small 
compared  to  kT  and  the  integrated  contribution  to  the 
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partition  function  is  well  represented  by  the  WKB  phase  shift. 


Fig.  7.  Plot  of  £=0  phase  shift  for  Li^CX^j;^).  Levinson’s  theorom 
.  c  g 

requires  a  threshold  value  q/tt  =  42.  The  rapid  fall-off, 
i.e.,  .3n/3e  «  0,  will  yield  a  large  negative  contribution  to 
the  molecular  partition  function.  This  curve  was  calculated 
using  the  WKB  theory  and  is  indistinguishable  from  the  exact 
quantal  results  (Shown  at  several  points)  except  at  energies 
within  a  few  cm“^  of  threshold.  Such  discrepancies  will  have 
a  negligible  effect  on  high  temperatures  equilibrium 
properties. 


Fig.  8.  Using  the  classical  theory,  we  show  the  contributions  to  q, 

•  f 

Li2(X^Eg)  at  6000°K.  These 
components  must  be  weighted  by  the  rotational  degeneracy  (2  +1 
to  obtain  the  total  internal  partition  function  Qi  = 

Z 

^UZ’ 

- total, - (B), - (M), - (F). 

Both  the  metastable  (M)  and  bound  (B)  components  are  always 
positive.  However,  the  free-continuum  (F)  phase  space  often 
produces  a  negative  contribution,  especially  for  high 
temperatures  and  for  shallow  potentials.  Only  alternate 
values  contribute  to  when  the  dimer  is  homonuclear. 


Fig.  9.  Contributions  to  q^  ^  for  Li2(b  at  6000°K.  See  Fig.  8 
for  explanation. 
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Fig.  10.  Phase  shift  ^qq  h  for  Jl=100  state  of  Li2(X^Z^).  Quantum 
theory  requires  a  threshold  value  ^qq  h  =  11  due  to  the  11 
true  bound  states  shown  in  Fig.  4.  The  rapid  rise  ton^  lOO^’^ 
18  at  e=1000  cm"^  is  a  result  of  the  7  metastable  states 
predicted  in  Fig.  4.  The  exact  quantum  results  (dotted)  show 
almost  step-like  increases  in  ^qq  by  modular-ir  at  each 
resonance  energy.  The  WKB  phase  shift  (solid  curve)  passes 
smoothly  through  this  region  and  joins  exactly  with  the 
quantum  phase  shift  above  the  rotational  barrier. 

Fig.  11.  Phase-shift  ^3Q/if for 2-=130  state  for  Li2(X^Zg).  Fig.  2 
shows  that  v-j  does  not  support  any  true  bound  states, 

and  the  threshold  value  must  be  ^3q  =  0.  A  rapid  positive 
increase  occurs  between  600  cm"^  and  2000  cm"^  where  the 
potential  supports  the  9  metastable  states  located  by  the 
step-like  structure  in  the  exact  quantum  phase  shift  (dotted). 

Fig,  12.  Phase  shift  n-j  i3q/^  for  S.=180  state  of  Li2(X^Zg). 

The  potential  v^  ^  indicates  that  neither 

bound,  nor  metastable  states  exist,  and 
Between  threshold  and  e=7000  cm”^  there  is  a  positive 
gradient,  9n/9c  >  0,  and  the  resultant  partition  function 

^1  180  “  '^l  180^^^  positive  for  temperatures  up  to 

10,000°K  and  higher.  Eventually  the  negative  gradient 
portions  of  ^gQ  will  tend  to  diminish  q^  until,  for 


definite  temperature,  there  is  a  perfect  cancellation,  since 
y^”de3n/3n  =  0  as  required  by  Levinson's  theorom. 


Fig.  13. 


Calculated  equilibrium  constant  K  =  -<B>  for  Li,  and  Na,. 

eg  t  t 

The  solid  curve  is  the  result  using  the  complete  partition 

functions  for  both  the  singlet  and  triplet  states.  The  dashed 

3  + 

curve  shows  the  effect  of  neglecting  the  b  state.  This 
deviation  will  be  manifest  as  a  residual  second  virial 
correction  which  must  be  applied  to  the  alkali  monomer  if  one 
hopes  to  reproduce  the  exact  equation  of  state  for  the  vapor. 


3 

Fig.  14.  Calculated  equilibrium  constant  Kp  (  E^)  for  the  triplet  state 
of  Li2.  The  solid  curve  is  obtained  by  assuming  Qjqj  =  and 
is  meant  to  mimic  the  qualitative  behavior  we  would  expect  for 
a  shallow  van-der-waals  type  dimer.  The  bound  state 
contribution  Kp(B),  shown  by  the  dashed  curve,  remains 
position  for  all  T,  but  the  combined  quantity  Kp  =  Kp(B)  + 
Kp(C)  predicts  a  negative  concentration  of  dimers  for  T 
5500°K.  This  effect  can  be  properly  described  in  terms  of 
modified  activity  for  the  monomer. 
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Barriers  to  Internal  Rotation  in  Inorganic  Species 


Stanley  Abramowitz 


ABSTRACT 

Barriers  to  interal  rotation  have  been  determined  using  infrared  and 
Raman  spectroscopy  for  some  inorganic  species.  The  determination  of  these 
barriers  will  be  described  for  a  high  barrier  (BCJI2SH)  a  medium  barrier 
(PFg,  AsFg  and  VFg)  and  a  low  or  zero  energy  barrier  (B(CH2)3). 


INTRODUCTION 


The  determination  of  the  energy  barriers  to  internal  rotation  and  the 
vibrational  energy  levels  associated  with  these  barriers  are  very  important 
for  the  computation  of  thermodynamic  functions.  This  paper  will  illustrate 
the  experimental  determination  of  barriers  to  internal  rotation  for  the 
species  BCJI2SH  (D),  XFg  where  X  =  P,  As  or  V,  and  8(0112)3.  These  represent 
a  high,  medium  and  low  (zero)  barrier  respectively.  The  techniques  of 
vibrational  spectroscopy  were  utilized  to  determine  the  relevant  vibrational 
energy  levels.  In  the  case  of  8(0113)3  a  comparison  of  the  third  law 
entropy  with  an  experimentally  determined  entropy  was  made. 

8CJI2SH 

The  8Cil2SH  molecule  and  its  deuterium  analog  have  a  planar  0^  symmetry. 

This  symmetry  has  7a'  and  2a"  modes  (1).  (It  should  be  noted  that  another 

Cj  structure  is  possible  with  a  plane  bisecting  the  CJ18CJI  angle. 

This  structure  has  6a'  and  3a"  vibrations.)’  The  vibrational  assignment  for 

the  8CJI2SH  and  8C1I2SD  species  is  given  in  Table  1.  (Where  8  and/or  CZ 

11  35 

isotope  effects  have  been  observed  the  frequency  given  is  for  the  8  C£2SH 
(D)  species.)  This  assignment  was  arrived  at  through  interpretation  of  the 
infrared  and  Raman  spectra  and  the  boron,  chlorine  and  hydrogen  vibrational 
isotopic  shifts.  The  vibrational  assignment  for  80)^2511  and  80^250  was  used 
to  calculate  the  force  field  for  these  species.  The  computed  force  field 
successfully  predicts  the  isotopic  splittings  for  these  species  (1).  The 
off  diagonal  force  constants  were  taken  from  8013  (2).  Other  off-diagonal 
force  constants  were  used  only  if  they  improved  the  overall  frequency  shift. 
The  syiiimetry  coordinates  and  force  field  and  the  assumed  geometric  structure 
for  the  8CJI2SH  (0)  species  can  be  found  in  reference  1. 


The  SH  (D)  torsional  mode  shows  a  hot  band  structure  in  the  infrared 
spectrum  with  0-1,  1-2,  2-3,  and  3-4  modes  observed  for  both  species.  Spectra 
observed  at  temperatures  higher  than  ambient  verified  the  assignments. 

These  modes  could  be  fit  with  either  of  two  potentials 

V  =  A(Z^  -  BZ^),  where  i  =  2(P  or, 

V  =  (V2/2)  (I-COS24.) 

This  is  because  the  barrier  is  quite  high  and  the  first  potential 

is  simply  (l-cosn(|))  with  a  truncated  power  series  substituted  for  cosnc)). 

4  2 

The  A(Z  -  BZ  )  which  is  most  useful  for  ring  puckering  modes  does  not 
require  a  knowledge  of  the  molecular  dimensions.  This  is  an  advantage 
for  BCJi2SH  (D)  since  experimentally  determined  bond  lengths  and  angles 
are  not  available.  The  V2/2  (l-cos24>)  which  requires  a  knowledge  of 
the  molecular  dimensions  is  a  better  model  for  the  molecular  motion  involved 
in  the  SH(D)  torsional  mode,  Vy. 

For  this  potential  the  moments  of  inertia  are  computed  from  the  assumed 
moleculer  dimensions.  Table  2  presents  a  summary  of  the  calculated  and 
observed  vibrational  frequencies  for  the  torsional  vibrational  modes  about 
the  B-S  bond.  As  one  can  see  the  barrier  is  quite  high.  The  more  physically 
meaningful  potential,  V=(V2/2)  (l-cos2<}>)  yields  a  barrier  of  about  4200  cm  ^ 
for  these  species. 

XFg  where  X=P,  As  and  V 

All  of  these  species  have  been  shown  to  be  trigonal  bipyramids  having 
structure.  The  F  atom  interchange  in  these  species  can  be  accomplished 
through  the  motions  along  the  E'  coordinate  which  can  lead  to  axial  equatorial 
interchange  through  a  intermediate  (3,4).  It  was  therefore  very  interesting 
to  have  obtained  the  Raman  gas  phase  spectra  of  these  species  in  order  to 
directly  observe  the  hot  band  structure.  This  vibrational  hot  band  structure 
must  be  interpreted  using  a  doubly  degenerate  harmonic  oscillator  function  (3,4). 
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The  Hamiltonian  for  axial -equatorial  interchange  is  given  by 
H  =  (a1i^/m)^/^ 

where  Hp  =  [1/2  (Px^  +  Py^)^  +  l/2{xV)-B(x^-3xy^)+C(x V)^] 

where  B  and  C  contain  a  term  which  is  a  scale  factor  involving  the  reduced 
mass  involved  in  the  vy  vibration  and  the  harmonic  force  constant  describing 
the  Vj  vibration.  One  can  produce  a  table  of  dimensionless  eigenvalues  for 
this  potential  (4),  where  the  barrier  for  F  atom  interchange  is  -(1/32)B  .  The 
important  feature  of  this  potential  is  that  it  requires  (n  +  1)  vibrational 
levels  where  n  is  the  vibrational  quantum  number. 

The  advantages  of  working  with  the  above  reduced  Hamiltonian  are  (1)  one 
need  not  calculate  an  explicit  reduced  mass  and  (2)  the^'e  is  one  less 
potential  constant  to  vary  than  in  the  explicit  form.  Thus  a  fit  of  the 
reduced  Hamiltonian  Hp  to  observed  energy  levels  will  yield  B,  C  and  the 
scale  factor  and  the  barrier  to  axial  -  equatorial  interchange. 

The  observed  and  calculated  transitions  for  PFg,  AsF^,  and  VFg  are 
given  in  Tables  3,  4,  and  5  respectively.  The  barriers  decrease  monotonically 
from  about  1050  to  475  cm"\  More  features  are  observed  in  VFg  then  AsF^  and 
more  in  AsF^  then  PF^  because  of  the  Boltzmann  factor.  The  vibrational 
frequency  decreases  from  175  to  133  to  108  cm  ^  in  going  from  PF^  to  AsFg  and 
VFg.  Therefore,  the  largest  number  of  vibrational  levels  are  populated  in 
VFg;  with  less  in  AsF^  and  even  less  in  PFg.  Further  features  of  the  dynamics 
of  axial  -  equatorial  fluorine  atom  interchange  can  be  found  in  the  references 
cited  (4). 

The  third  law  entropy  of  B(CH2)2  computed  from  the  available  vibrational 
assignment  and  inertial  parameters  assuming  free  rotation  of  the  CH^  groups 
is  not  in  agreement  with  the  experimentally  measured  entropy  (6,7).  (The 


-3- 


difference  is  such  that  any  barrier  to  rotation  would  increase  the  discrepancy). 
Since  the  Raman  gas  phase  spectrum  for  this  species  had  not  been  previously 
observed  a  reinvestigation  of  the  vibrational  spectrum  was  initiated  in 
order  to  better  characte'^ize  the  vibrational  assignment.  Boron  trimethyl 
has  a  planar  BC3  skeleton.  Conventional  point  group  theory  would  indicate 
an  average  symmetry  for  freely  rotating  CH^  groups  with  20  normal  modes. 
However,  the  observed  infrared  and  Raman  gas  phase  spectra  showed  fewer 
features.  The  observed  spectra  indicated  that  the  6^24  symmetry  group  based 
upon  the  molecular  symmetry  (MS)  group  theory  first  proposed  by  Lonquet 
-  Higgins  (8,9)  was  necessary.  Since  the  discussion  of  the  vibrational 
structure  in  terms  of  the  6224  group  and  its  relationship  to  has  not 
been  previously  given  a  brief  discussion  of  the  vibrational  analysis  follows. 
Finally  a  program  developed  to  compute  the  thermodynamic  functions  for 
anharmonic  oscillators  was  utilized.  Somewhat  better  agreement  between 
the  computed  third  law  entropy  and  that  experimentally  determined  is 
obtained. 

The  vibrational  modes  of  B(CH2)2  shown  in  Table  6.  For  a 
symmetry  with  freely  rotating  methyl  groups  one  expects  20  normal  modes. 

The  Aj ,  E'  and  species  are  Raman  active  while  the  and  E'  species  are 
infrared  active.  If  one  considers  the  0^24  molecular  symmetry  group 
there  are  15  normal  modes  of  which  the  A-j,  I,  G  and  E2  are  Raman  active  while 
the  A^,  I,  G  and  E2  are  infrared  active.  Table  5  shows  a  correlation  diagram 
between  and  G224  symmetries.  It  should  bo  noted  that  6^24  requires 
8  Raman  active  modes  and  nine  infrared  active  modes  while  symmetry 
requires  13  Raman  active  modes  and  11  infrared  active  modes.  The 
significant  feature  of  the  MS  group  is  tl\e  relative  simplicity  of  the 
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expected  vibrational  spectrum.  The  observed  Infrared  and  Raman  vibrational 
spectra  confirm  the  MS  symmetry  group  as  the  better  representation  of  this 
molecule's  symmetry.  The  BC^  asymmetric  stretching  frequency  for  or 
for  6224* symmetries  respectively  shows  a  isotope 

shift  of  about  30  cm"^.  This  shift  Is  to  be  expected  for  a  planar  BC^  skeleton. 
(The  corresponding  mode  in  the  boron  trihalides  also  show  a  similarly  large 
Isotope  shift.)  Three  out  of  the  four  modes  have  the  same  P-R  separation 
which  is  to  be  expected  for  these  parallel  transitions.  As  Is  to  be  expected 
these  modes  are  not  in  the  Raman  spectrum.  The  three  A^  modes  at  2916,  1291, 
and  678  cm~^  are  all  polarized  and  are  absent  in  the  Infrared  spectrum  as 
required  by  the  selection  rules.  The  two  6  fold  degenerate  vibrations  and 
v^2  have  rather  broad  unstructured  contours  in  both  the  Infrared  and  Raman 
spectra.  In  summary  the  Infrared  and  Raman  spectra  are  consistent  with  a 
6^24  molecular  symmetry  assignment.  There  are  not  enough  vibrational  bands 
observed  to  justify  a  assignment  for  the  B(CH2)2  molecule.  A  summary  of 
the  vibrational  assignment  Is  given  In  Table  7. 

A  calculation  of  the  third  law  entropy  of  8(0113)3  can  be  made 
using  the  rigid  rotor  harmonic  oscillator  approximation  and  by  assuming 
free  rotation  of  the  three  methyl  groups  about  the  B-C  bonds.  The 
vibrational  assignment  given  in  Table  7  Is  used  together  with  the  bond 
distances  and  angles  determined  using  electron  diffraction  techniques  by 
Bartell  and  Carroll  (10).  An  entropy  of  67.83  cal/K  at  199.92  K  Is  computed 
assuming  free  rotation  of  the  CH^  groups.  This  compares  with  the  measured 
value  of  68.29  +  0.10  cal/K  (1).  The  fact  that  the  measured  entropy  is 
greater  than  that  computed  using  spectroscopic  data  Indicates  essentially 
free  rotation  of  the  methyl  groups.  Any  significant  barrier  to  the  rotation 
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of  the  methyl  groups  results  in  a  lowering  of  the  calculated  entropy  and 
therefore  an  increase  in  the  discrepancy  between  the  measured  and  computed 
entropy. 

A  computer  program  has  been  written  to  compute  the  effect  of 

anharmonicity  on  the  vibrational  contribution  to  the  entropy.  In  this 

computation  all  vibrational  states  were  assumed  to  have  a  1%  anharmonicity 

(other  values  of  the  anharmonicity  can  be  chosen  for  each  vibration).  This 

computer  code  obtains  a  state  sum  by  counting  those  levels  whose  partition 

"8  "ll 

coefficient  exceeds  a  chosen  value,  in  this  case  10  .  (A  choice  of  10  does  not 
effect  the  computed  entropy.)  An  increase  of  0.03  e.u.  at  199.92K  is  obtained  in 
this  manner.  The  projected  difference  increases  for  higher  temperatures. 

This  program  will  be  especially  useful  in  computing  estimated 
errors  in  entropy  caused  by  vibrational  anharmonicity.  It  will  therefore 
be  useful  for  comparing  third  law  (statistically  calculated)  entropies 
with  those  derived  from  thermodynamic  measurements.  By  identifying 
expected  differences  the  evaluation  of  thermodynamic  data  particularly 
for  high  temperature  species  should  be  expedited. 
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Table  1.  Vibrational  Assignment  BCtg  SH  and  BCtg  SD 


Vibrational 

Assignment 

Approximate 

Description 

BCt2SH 

6Ct2SD 

'’I 

SH  (D)  str. 

2602 

1863 

''2 

BC)i2  asymm.  str. 

994 

999 

^3 

BS  str. 

900 

900 

A' 

SH(D)  bend 

779 

599 

''5 

BC«,2sym  str. 

461.8 

459 

CtBCt  bend 

255 

255 

'’7 

CtBS  bend 

230 

227 

00 

■> 

< 

Ct  0.0. p. 

474 

475 

'^9 

SH(D)torsion* 

386 

283 

★ 

The  0-1  frequency  is  given. 

Hot  band  transitions  are  given  in  Table  2. 


j 

i 
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Table  3.  Observed  and  Calculated  Transitions  for  PFg 


Transitions 

Symmetry 

Calculated 

(cm"^) 

B  =  0.07 

Calculated 

(cm'^) 

B  =  0.075 

Observed' 

(cm’b 

Raman 

0  -*■  1 

Aj  -V  E' 

175.0 

175.0 

175 

1  2 

H  •*  '‘1 

161.8 

159.8 

159 

1  2 

E'  ->  E’ 

177.9 

178.0 

178 

2-^3 

A]  -►  E’ 

164.8 

163.5 

164 

2-^3 

Aj  *Ai 

192.7 

194.0 

194 

2  ->  3 

Aj->Aj 

198.9 

199.4 

b 

2  3 

E'  -»■  E‘ 

148.8 

145.3 

b 

2-^3 

Aj.Aj 

177.5 

175.7 

c 

2  ->  3 

Ui 

182.9 

184.3 

180 

Barrier 

height 

(cm"^) 

1139 

995 

®  The  additional  peaks  at  168  cm"^  and  187  cm”^  probably  correspond  to  OP  and 
RS  maximum,  respectively. 

^  Not  observed. 


c 


This  feature  is  probably  obscured  by  0  -*■  1  transition 


Table  4.  Observed  and  Calculated  vy  Transitions  for  AsFj. 


Transition 

Symmetry 

Calculated 

(cm"’) 

B  =  0.07 

Calculated 

(cm"’) 

B  =  0.075 

Observed 

(cm"’) 

Raman 

0  H.  1 

A|  -  £■ 

132.8 

132.8 

133 

1  X.  2 

E*  -  Aj 

122.8 

121.2 

121 

1  2 

E‘  -  E' 

135.0 

135.1 

136 

<o 

f 

CM 

A|  -  E* 

125.1 

124.1 

125 

2  X.  3 

Aj-AJ 

146.2 

147.2 

148 

2-3 

Aj-Aj 

151.0 

153.7 

157 

2-3 

E'  -  E' 

112.9 

110.2 

113 

2-3 

E‘  -  Aj 

134.7 

133.3 

a 

2-3 

E'  -  Aj 

138.8 

139.9 

141 

3-4 

A^-A- 

117.6 

116.2 

118 

3-4 

E‘  -  E' 

126.7 

125.3 

129 

3-4 

E'  -  E’ 

161.7 

165.6 

166? 

3-4 

Aj-A^ 

96.5 

93.1 

b 

3-4 

A^  -  E* 

’  105.4 

102.2 

b 

3-4 

A|  -  E’ 

140.5 

142.5 

143 

3-4 

Aj-A| 

91.7 

86.6 

b 

3-4 

A|  -  E' 

100.7 

95.7 

b 

3-4 

A]  -  E’ 

135.8 

135.9 

c 

Barrier 

height 

(cm"’) 

... 

864 

755 

*  Obscured  by  0  1. 

^  Not  observed. 

'  Obscured  by  1  ►  2(E'  ->  E’). 


Table  5.  Observed  and  Calculated  Transitions  for  VFg. 


Transition 

Symmetry 

Calculated 

(cm~^) 

B  =  0.08 

Calculated 

(cm'b 

B  =  0.09 

Observed  (err 
Infrared® 

Raman^ 

0  1 

A'  - 

E' 

107.4 

197.3 

— 

108 

1  ^  2 

E'  - 

A' 

''l 

96.8 

94.0 

91.9 

93 

1  ^  2 

E'  - 

E' 

109.4 

109.4 

109.4 

108 

2-^3 

A|- 

E’ 

99.3 

97.0 

98.1 

100 

2  3  R 

Aj- 

A' 

''l 

119.5 

119.5 

c 

120 

2  3  R 

Ai- 

A' 

''l 

126.6 

131.6 

c 

126,120 

2^3 

E'  - 

E' 

86.7 

81.6 

85.8 

5,83 

2-3 

Aj  - 

A' 

''l 

107.0 

104.2 

103.6 

— 

2-3 

E'  - 

A' 

114.0 

116.2 

— 

114 

Barrier 

height 

(cm"^ ) 

— 

539 

428 

— 

®  Reference  5. 

^  There  are  additional  Raman  features  which  can  be  attributed  to  3  -»■  4 
transitions.  Limitations  in  our  calculations  did  not  enable  us  to 
determine  these  transitions. 

^  Only  Raman  active. 


Table  6.  Correlation  Table  for  6^24  and  for  B(CH2)2 


Point  Group 


^3h 


Molecular  Symmetry  Group 


'324 


Species 

Vibration 

Activity 

Mode 

^1 

1 

R(P) 

CH3  stretch 

2 

CH3  deformation 

3 

BC3  stretch 

a; 

4 

lA 

CH3  stretch 

5 

CH3  deformation 

6 

CH3  rock 

II 

^2 

7 

IR 

CH3  stretch 

8 

CH3  deformation 

9 

CH3  rock 

10 

BC3  deformation 

E* 

11 

IR,R 

CH3  stretch - 

12 

CH3  stretch 

13 

CH3  deformation— - 

14 

CH3  deformation^  ' 

15 

16 

BC3  stretch.. 

CH3  rock 

17 

BC3  deformatiori  . 

E" 

18 

IR.R 

CH3  stretch  j 

19 

CH3  deformation 

20 

CH3  rock 

Vibration  Species  Activity 


1 

2 

3 


4 

5 

6 


7 

8 
9 

10 


13 


s.jjv-'  I*' 


R(P) 


A4  lA 


A3  IR 


IR,R 


IR,R 


Eg  IR.R 


Table  7. 

Vibrational  Assignment  of  B(CH3)3 

According  to  633^  MS  Group 

Species 

Vibration 

JR 

Raman 

Inactive  (Estimate! 

1 

2916 

2 

1291 

3 

678 

4 

2915 

5 

1400 

6 

900 

h 

7 

2914 

8 

1303 

9 

975 

10 

326 

I 

n 

2984 

2985 

12 

1450 

1445 

G 

13 

867 

868 

^2 

14 

1184® 

1162 

1154. 

15 

310 

314 

®  isotopic  splitting 


Studies  of  the  Reactions  of  M(g)  +  SO2 
S.  Abramowitz  and  L.U.  Sleek 

The  matrix  isolated  infrared  spectra  of  the  products  of  reactions  of 
M(g)  where  M=  ^Li,  ®Li,  Na  and  K  with  S02(g)  have  been  observed.  Most  of 
the  effort  has  been  directed  towards  the  Li  +  SO2  system.  Absorption  bands 
in  these  spectra  have  been  assigned  to  a  symmetric  and  asymnetric  stretching 
mode  of  SO2".  There  are  other  absorption  bands  in  the  spectra  which  are  due 
to  M-SO2  motions.  Isotopic  studies  utilizing  S^^02  and  S^^O^^O  as  well  as  ^Li 
and  ^Li  are  expected  to  elucidate  the  structure  of  the  MSO2  species.  In  particular 
we  are  hopeful  that  our  results  can  indicate  whether  the  M  is  bonded  to  the 
S  or  O2. 

The  results  of  this  study  will  enable  the  calculation  of  the  thermodynamic 
properties  of  the  MSO2  species.  Species  such  as  these  may  be  important  in 
the  corrosion  of  turbine  blades  by  high  temperature  gases  containing  sulfur 
oxide  species  and  trace  metals. 


